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Abstract. The motion of an elastic solid inside of an incompressible viscous fluid is 
ubiquitous in nature. Mathematically, such motion is described by a PDE system that 
couples the parabolic and hyperbolic phases, the latter inducing a loss of regularity which 
has left the basic question of existence open until now. 

In this paper, we prove the existence and uniqueness of such motions (locally in time), 
when the elastic solid is the linear Kirchhoff elastic material. The solution is found using 
a topological fixed-point theorem that requires the analysis of a linear problem consisting 
of the coupling between the time-dependent Navier-Stokes equations set in Lagrangian 
variables and the linear equations of elastodynamics, for which we prove the existence of a 
unique weak solution. We then establish the regularity of the weak solution; this regularity 
is obtained in function spaces that scale in a hyperbolic fashion in both the fluid and solid 
phases. Our functional framework is optimal, and provides the a priori estimates necessary 
for us to employ our fixed-point procedure. 

1. Introduction 

We are concerned with establishing the existence (and uniqueness) of solutions for the 
equations of motion of linearly elastic solids moving and interacting with an incompressible 
viscous fluid, with the natural conditions of continuity of the velocity fields and normal 
components of the stress tensors along the moving interface between the two materials. 

The analysis of interacting fluid-structure problems has been the subject of active research 
since the late nineties. As of now, only the question of the possible motion of a solid inside 
of a viscous flow, in which the solid is either rigid or consists of a finite number of modes, 
has been settled. In [T^J], existence and uniqueness (locally in time) of smooth solutions 
has been obtained using a Lagrangian framework, for the rigid body case, provided that 
the rigid disk is sufficiently heavy. In ^U], for the same problem, but with an arbitrary 
number of rigid solid bodies, existence of at least one weak solution has been established 
in an Eulerian formulation by a global variational approach; their result holds for all time 
in two space dimensions as long as no collisions occur between solids or with the boundary, 
and is local in time for the three dimensional case. In by generalizing the methods 
of jlUj . the case of an elastic body following the linear Kirchhoff law, with the important 
restrictions of allowing only a finite number of modes, and a relaxation of the continuity 
of the normal stress along the boundary of the solids, has been considered. The above list 
of references for contributions to this area is by no means exhaustive; see for instance [7], 
|14j . |17j . Note also that the related problem of the free fall of a rigid body in a Stokes 
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flow in the full space has been considered in [22], for the stationary case, and in |19j for the 
stationary as well as the time-dependent case. 

More recently, the interaction of a viscous incompressible flow with an elastic plate (with- 
out the restriction of a finite number of modes), whose constitutive law comprises a parabolic 
hyperviscosity term in the plate, has been studied in [5]. We remark that this additional 
hyperviscosity term is of crucial importance in that study. (Note also that two dimensional 
plate models that approximate thin three dimensional structures usually contain fourth- 
order operators arising from bending stresses, whereas models of elastic solids have only 
second-order operators; as such, plate models can provide better a priori control for the 
motion of the material interface.) 

In the steady-state situation in which both phases are governed by elliptic operators, JB] 
has obtained an existence result (for the case that the solid follows the nonlinear Saint 
Venant-Kirchhoff law) by the use of a fixed-point method that iterates between fluid and 
solid phases. This approach is indeed natural for the steady-state problem since the analysis 
can make use of elliptic regularity theory. For the dynamic problem, however, such an 
iteration procedure appears to fails because of a consequent loss of regularity induced by 
either a fluid-solid-fluid iteration or a solid-fluid-solid iteration. This loss of regularity is 
due to the fact that hyperbolic and parabolic systems do not have the same regularity 
requirements and properties, which is in fact the heart of the difficulty in the coupling of 
the two phases. 

Whereas the coupling between the Navier-Stokes equations and the linear Kirchhoff law 
is perhaps the most fundamental problem to consider in regards to the motion and inter- 
action of an elastic body in a viscous incompressible fluid, none of the methods that have 
been developed to date can handle this system, mostly because of the differences between 
parabolic and hyperbolic regularity, i.e. in both the requirements on the function spaces 
for the prescribed data, as well as the functional framework of the solution space. 

We now come to the formulation of the problem. The motion of the fluid is described 
by the time-dependent incompressible Navier-Stokes equations, while the deformation of 
the solid body is governed by the linear Kirchhoff equations. The two models are coupled 
along the moving material interface by imposing the continuity of the normal component 
of the stress tensors as well as the particle displacement fields. From the point-of-view of 
mathematical analysis, the Navier-Stokes equations are traditionally studied in the Eulerian 
(or spatial) description, while the elastic body is studied in the Lagrangian (or material) 
frame. Because the material interface is fixed in the Lagrangian representation, we shall 
study this problem entirely in material coordinates. This Lagrangian framework also has 
the advantage of keeping the hyperbolic problem (where the loss of regularity occurs) linear, 
which is of paramount importance here. Note, however, that a semi-linear elastic system, 
as for some plate or shell models (see for instance 0), can be handled without any difficulty 
by our methodology. The question of existence for the case of a quasilinear elasticity law 
can also be obtained (and shall be addressed in later work), requiring a smoother functional 
framework leading to more compatibility conditions at the origin. 

Let us now set the equations. Let C R 3 denote an open, bounded, connected and 
smooth domain with smooth boundary d£l which represents the fluid container in which 
both the solid and fluid move. Let Q s (t) C Q denote the closure of an open and bounded 
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subset representing the solid body at each instant of time t G [0, T] with W (t) := Vl/Vl s (t) 
denoting the fluid domain at each t € [0, T\. Note that in our analysis VI s (t) is not necessarily 
connected, which allows us to handle the case of several elastic bodies moving in the fluid. 

Remark 1. If a function u is defined on all of VI, we will denote u$ = u 1^/ and u s = u 1q« . 
This allows us to indicate from which phase the traces on 



r(o) := vif(o) n vi s {o) 

of various discontinuous terms arise, and also to specify functions that are associated with 
the fluid and solid phases. 

For each t 6 (0, T], we wish to find the location of these domains inside Vl, the divergence- 
free velocity field vJ (t, •) of the fluid, the fluid pressure function p(t, •) on VI* (t), the fluid 
Lagrangian volume-preserving configuration r]f(t,-) : f2^(0) = VIq — ► Vl f (t), and the elastic 
Lagrangian configuration field r] s (t, •) : VI s (0) = VIq —>■ VI s (t) such that 

Vl = v s (t,%)U V f(t,Vl f ), (1.1a) 

where 

rj{{t,x) = U f (t,r l f(t,x)), (1.1b) 

and solves the Navier-Stokes equations in VI? (t) 

u{ + (u f -V)u f = divT f + f f , (1.1c) 

drvV = 0, (l.ld) 

with 

T f = u Defu f -pi, (Lie) 
and rj s solves the elasticity equations on VI s (0) 

f = d\vT s + f s , (l.lf) 

with 

T s = A Trace(V?? s - /)/ + ft (Vr/ S + V V sT -21), (l.lg) 

and where the equations are coupled together by the continuity of the normal component 
of stress along the material interface T(t) := Vl s (t) D Vl f (t) expressed in the Lagrangian 
representation on 

T := r(0) 

as 

T s N = [T f o n 1 ] [(Wy 1 N] , (l.lh) 
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and the continuity of particle displacement fields along Tq 

r] f = V s , (1-li) 
together with the initial conditions 

u(0,x) = u {x) , 
77(0, x) = x , 

and the Dirichlet (no-slip) condition on the boundary dQ of the container 

u f = 0, (1.11) 

where v > is the kinematic viscosity of the fluid, A > and /i > denote the Lame 
constants of the elastic material, N is the outward unit normal to To and Def u is twice the 
rate of deformation tensor of u, given in coordinates by u l ,j +u J ,j. All Latin indices run 
through 1,2,3, the Einstein summation convention is employed, and indices after commas 
denote partial derivatives. 

We now briefly outline the proof. As the solid and fluid phases are naturally expressed 
in the Lagrangian and Eulerian framework, respectively, we begin by transforming the 
fluid phase into Lagrangian coordinates, leading us to the system of equations (|3.2|) of 
Section |31 This system of PDE is both parabolic (in the fluid) and hyperbolic (in the 
solid) in character; hence, one of the fundamental difficulties that must be overcome is 
an appropriate functional framework accommodating both features. Sections and are 
devoted to the setting of our functional framework, which appears to be of hyperbolic-type 
in both solid and fluid phases, and is necessitated by the estimate of the elastic energy. 
This hyperbolic scaling in turn requires the initial data to possess more regularity, and 
thus produces more compatibility conditions in the fluid phase than if a parabolic scaling 
were used (as seen in the statement of the existence theorem in Section 0). Whereas the 
choice of working in Eulerian or Lagrangian variables may seem arbitrary, at the level of the 
functional framework, it appears that the problem truly requires this hyperbolic functional 
framework for both phases, regardless of the choice of spatial or material coordinates. 

In order to solve (|3.2|) . we use a fixed-point approach, where we solve the linear system 
([71]) for the La grangian velocity w, the coefficients a^rj) coming from the flow map rj of 
a given velocity v. The study of the regularity of the solutions to this problem, which 
constitutes the main part of this paper, is given in Sections H and CHI It appears that 
the regularity theory for ()7.1|) cannot be obtained directly by solving the problem with the 
actual coefficients oc(^). In Section|H]we explain the smoothing process for the problem: we 
introduce smoothed velocity fields v n which provide us with smoothed coefficients a{(rj n ) 
(which we denote generically by v and a). We also present two versions of what we term 
the Lagrange multiplier lemma (which associates a pressure function to the weak solution) 
that will be of basic use throughout this paper. 

We study in Section ED the existence of weak solutions w to (|7,lj) (with regularized co- 
efficients), as the limit of penalized problems. Whereas these penalized problems are not 
necessary merely to obtain existence of weak solutions, they are of paramount importance 
in getting the appropriate regularity results for wt and wu, the primary reason being that 
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the pressure associated to (|7.1|) with the Dirichlet boundary condition cannot be obtained 
simply from the variational form of the problem, and requires the study of the time dif- 
ferentiated problem in order to get more information on wt (which would need to be in 
L 2 (0, T; H^ 1 (O; M 3 )) for the Lagrange multiplier lemma). Unfortunately, this time differ- 
entiated problem contains p o rj in its formulation, which leads to a circular argument, and 
thus explains the need for the penalized problem. We then obtain the regularity for the 
problem by the energy inequality for wu and some difference-quotient inequalities for wt 
and w carried-out in Lagrangian variables in a neighborhood of the interface Tq. This, 
in turn, provides us with an estimate for the trace of w and wt on To, which after a re- 
turn to the Eulerian variables for the fluid phase, immediately provides the regularity in the 
fluid domain. The regularity in the solid phase is then obtained in a straightforward manner 
from elliptic regularity and the already-obtained trace estimate. We note that the estimates 
proved at this stage blow-up as the regularized coefficients tend to the true coefficients, i.e., 
as the regularization parameter tends to zero. 

For this reason, in Section 1101 we obtain a different set of estimates (founded upon 
interpolation inequalities) for the solutions of the regularized problems, and conclude that 
the norms of the regularized solutions are actually uniformly bounded in the appropriate 
spaces, which thus provides by weak convergence, a solution to (|7.1|) with the appropriate 
a priori estimates. 

Finally, we conclude the proof of the existence theorem in Sections ^2 and El by means of 
the Tychonoff fixed-point theorem. Although it might be possible to employ the Schauder 
theorem instead, it appears that the strong convergence requirements of the Schauder the- 
orem are not very convenient to write and are, in particular, unnecessary for the use of the 
Tychonoff theorem. 

Uniqueness is obtained with more regularity on the forcing functions and initial data in 
Section 1131 in order to get information on the second-time derivative of the pressure qu of 
the solution built, that we do not have with the assumptions of Theorem 15. II 



2. NOTATIONAL SIMPLIFICATION 

Although a fluid with a Neumann (free-slip) boundary condition indeed obeys the con- 
stitutive law (jl.lej) . it turns out that the notation is substantially simplified (particularly 
in Section E3 wherein we analyze the twice differentiated-in-time problem in Lagrangian 
coordinates) if we replace (jl.lej) with 

T f = uVu f -pl; (2.1) 

this amounts to replacing the energy / Def-u-^ : Def v by / Vu-^ : Vi>, which is an 

equivalent form when = on dQ. due to the well-known Korn inequality. Henceforth, we 
shall take (|2.1j) as the fluid constitutive law. 



3. Lagrangian formulation of the problem 

In regards to the forcing functions, we shall use the convention of denoting both the fluid 
forcing ff and the solid forcing f s by the same letter /. Since // has to be defined in fi 
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(because of the composition with 77), and f s must be defined in Qq, we will assume that the 
forcing / is defined over the entire domain f2. 
Let 

a(x) = [Vr/(£)r\ (3.1) 
where (Vrjf(x))j = d(r}* Y/dx^(x) denotes the matrix of partial derivatives of 7/*. Clearly, 
the matrix a depends on 77 and we shall sometimes use the notation o*-(r/) to denote the 
formula l|3.1|) . 

Let v = uot] denote the Lagrangian or material velocity field, q = por] is the Lagrangian 
pressure function (in the fluid), and F = f* o r/f is the fluid forcing function in the material 
frame. Then, as long as no collisions occur between the solids (if there are initially more 
than one) or between a solid and d£l, the system (|l.lj) can be reformulated as 

7] t = v in (0,T) x O, (3.2a) 

v\ - u(aja k v l , k ), 3 +(aj?), fc = F i in (0,T) x fij , (3.2b) 

afV,*, = in (0,T) x fij , (3.2c) 

v\ - [c ijkl f v k 4 l, = f in (0, T) x % , (3.2d) 

j 

v v\ k a k l a) l Nj-qa\Nj = t" kl I v k ,1 Nj on (0,T)xr , (3.2e) 

^VjgflJfU;! 3 ) a.e. in (0,T), (3.2f) 
1; = no on O, x {t = 0} , (3-2g) 

77 = Id on J! x{t = 0), (3.2h) 

where N denotes the outward-pointing unit normal to Tq (pointing into the solid phase), 
and 

c ijkl = X5 ij 5 kl + ^gikgjl + fiUfiih} _ 

Throughout the paper, all Greek indices run through 1, 2 and all Latin indices run through 
1,2,3. Note that the continuity of the velocity (jl.lij) along the interface is satisfied in the 
sense of traces on To by condition (|3.2fj) . whereas the continuity of the normal stress along 
the interface is represented by (|3.2ej) . 

Remark 2. The case in which the viscosity or Lame coefficients are variable functions 
depending on x £ fi and satisfying the usual assumptions, can be handled by our methodology 
without any supplementary mathematical difficulties. 

4. Notation and conventions 

We begin by specifying our notation for certain vector and matrix operations. 

We write the Euclidean inner-product between two vectors x and y as x ■ y, so that 
x • y = x 1 y % . 

The transpose of a matrix A will be denoted by A T , i.e., (A T Yj = A\. 

We write the product of a matrix A and a vector b as A b, i.e, (A b) 1 = A l -b> . 

The product of two matrices A and S will be denoted by A -S, i.e., (A-S)^ = A\ S k . 
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The trace of the product of two matrices A and S will be denoted by A : S, i.e., 
A : S = Tvace(A ■ S) = A) S{. 

For s > and a Hilbert space (X, \\ ■ \\x), H S (Q;R 3 ) denotes the Sobolev space of R 3 - 
valued functions with s distributional derivatives in L 2 (Q;R 3 ), while L 2 (0,T; X) denotes 
the equivalence class of functions which are measurable and have finite || • 11^2 -norm, where 

ll/ll! 2 (o,T;X) = /o r H/WII^- 

We also set F^(^;M 3 ) ={«£ H 1 (Q f ;R 3 )\ u = on dfl} . 
For T > 0, we set 

Vf{T) = {w G L 2 (0,T;H 2 (Q f ;R 3 )) | w t G L 2 (0, T; fl^fij; M 3 )) 

w« GL 2 (0,T;L 2 (^;M 3 ))}, 

V*(T) = {u; G L 2 (0,T;H 3 (Q f ;R 3 )) \ w t G L 2 (0, T; # 2 (^; M 3 )) | 

y s 2 (T) = {u> G L 2 (0,T;F 2 (^;M 3 )) | w t G L 2 (0, T; tf 1 ^; R 3 )) | 

Wii a 2 (0,T:L 2 (^;l 3 ))} , 
K 3 (T) = {w G L 2 (0,r;F 3 (^;M 3 )) | w t G L 2 (0, T; # 2 (5^; M 3 ))| 

«fc GL^O.T;]? 1 ^;! 3 ))}. 

We will solve (|3.2|) by a fixed point method, set in an appropriate subset of Vf (T) x 1/ 3 (T). 

We assume in what follows that v G V 3 (T) is given in such a way that the matrix a^if) 
associated to the flow rj of this velocity field v is well-defined. 
We then introduce the space (of weak solutions) 

V v ([0,T]) ={w G L 2 (0, T;L 2 (f2;M 3 )) | f w G L 2 (0, T; ff 1 (il;!R 3 )), 

JO 

if G L 2 (0, T; ^(JIq ; M 3 )), af = in [0, T] X w = 0ondQ}. 

Note that we impose the condition / w G L 2 (0, T; M 3 )) to ensure continuity of 

Jo 

the displacement field, in the sense of traces, between the solid and fluid phases along the 
interface To- We will also denote for t G [0, T] 

V v (t) = {V G H^(n;R 3 ) | = in o£} . 

Furthermore, we will need the space 

W([0, T]) ={w G L 2 (0,T;L 2 (S1;IR 3 )) | / iu G L 2 (0, T; ff^Jl; R 3 )), 

Jo 

w G L 2 (0, T; i7 1 (il f ; M 3 )) w = 0on 90}, 

with the "divergence-free" constraint removed. 

In order to specify the initial data for the weak formulation, we introduce the space 

L 2 divJ = {tpe L 2 (n;R 3 ) | div^ = 0inO£, ip • N = on dtt} , 
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which is endowed with the L 2 (Q;R 3 ) scalar product. 

The space of velocities, Xt, is defined as the following separable Hilbert space: 

X T = {ue L 2 (0,T;H^(n;R 3 )) \ {u^ J u s ) G Vf(T) x V 3 (T)} , (4.1) 

endowed with its natural Hilbert norm 

1 1 1 1 iiii2 1 1 1 1 2 1 1 1 1 2 

\\ U \\X T = \\ U \\L2(0,T;m(nm) + ll"ll i 2 (0)T . jff 3 (n /. K 3 )) + II U *IIl2(0,T;//2(^;R3)) 

l|/,|2 ||||2 ||i|2 

+ II / W llL 2 (0,T;tf 3 (ng;R 3 )) + \\ u \\ L 2 (0,T;tf 2 (ng;M 3 )) + ll n * IIl 2 (0,T;#^o; m3 )) 

J 
. II i|2 
+ H U ttllL2( 0)T;H l( n /. R 3^ ■ 

The existence of solutions to (|3.2|) will be obtained in the following separable Hilbert 
space 

Yt = {(u,p) eX T x L 2 (0,T;H 2 (nl;R))\ p t G L 2 (0, T; H 1 ^; R))} , 
endowed with its natural Hilbert norm 

ll(«»P)lly T = IMIx T + IIPllia (0) T;H3(n^)) + ll P *llL2(0,T;i?i(n^M)) ' 

Remark 3. Note well that our method does not require any a priori knowledge of the regu- 
larity of the second time derivative of the pressure function pui this is due to the Dirichlet 
boundary condition on dQ as well as the Lagrangian representation of the problem that we 
employ. 

We shall also need L°°-in-time control of certain norms of the velocity, which necessitates 
the use of the following closed subspace of Xt~- 

W T = {u£X T \ u tt G L°°(0,T;L 2 (ft;IR 3 )), u t G L°°(0, T; H 1 ^; R 3 )), 

u G L°°(0, T; H 2 (%; R 3 )), f u G L°°(0, T; H 3 (n s ; R 3 ))} , 

Jo 

endowed with the following norm 



n2 iiii2 II ||2 || / ||2 

\ u \\Wt = \\ U \\X T + INt|lL<x,( 0j T;L2(Q ; R3)) + || / u \\L°°(0,T;H3(n s ; 

J 



i || m2 _| || ||2 

+ H n llL°°(0,T;H 2 (Qg;]R 3 )) + II U * III, 00 (O.T;// 1 (Og 

For some of our estimates, we will also make use of the space 

Zt = {(u,p) G Wt x L 2 (0,T;If 2 (ft£;R))| p t G L 2 (0, T; fl^fiJjR))} , 
endowed with its natural norm 

II0,P)II! T = ||«||wt + \\P\\ 2 L 2 iQjT . H 2 {u f m) + IIPt& (0jTi Hi ( n/.R)) • 

Throughout the paper, we shall use C to denote a generic constant, which may possibly 
depend on the coefficients v, A, /i, or on the initial geometry given by Q and S7q (such as 
a Sobolev constant or an elliptic constant). Similarly, we will denote by C(M) a generic 
constant which depends on the same variables as C as well as on M (which is a variable 
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defined in the next section) and II u o||^5(q/. K 3)> ||/(0)||#3(n ; R3) an d the fixed time T for which 
the forcing functions are defined. We note that these constants do not blow-up whenever 
the quantities they depend upon remain finite. 

For the sake of notational convenience, we will also write u(t) for u(t, ■). 



5. The main theorem 

Theorem 5.1. Let O C M 3 be a bounded domain of class H 3 , and let Qq be an open set 
(with a finite number > 1 of connected components) of class H A such that fip c O. Let us 
denote n f = Q n (Hf) c . Let v > 0, A > 0, /x > be given. Let 

f G L 2 (0,f;H 2 (Q;R 3 )), f t G L 2 (0, T; F 1 ^; M 3 )), f tt G L 2 (0, f ; L 2 (Q; R 3 )), (5.1a) 

/(0)£ff 3 (!l;I 3 ). (5.1b) 

Assume that the initial data satisfies 

u g # 5 (q£ ; r 3 ) n H\n s o; r 3 ) n FoHfi; r3 ) n Lj ivJ 

as well as the compatibility conditions 

[Vuq ./V]tan = on Fq, w\ = on dQ, uAu^ — Vqo = on Tq, (5.2a) 

MVw{ NY + U[ U f ,l (O, fc of) t (0) ] ^)il]tan 

= [(<? iia 4>l N i)i=i]*m onT , (5.2b) 

with qo G H 4 (Qq;R) defined by 

A go = div/(0) + (a?) t (0)t4,i *n ^o> ( 5 - 3a ) 
go = z,[V?4 AT] • JV on T , (5.3b) 

|^ = /(0) • iV + i/Auo • N on 90, (5.3c) 

anrfwi G i^(0;K 3 ) niT 3 (Qg;K 3 ) niI 3 (n£;R 3 ) denned 6y 

wi = uAu - Vq + /(0) in Oq (5.4a) 
wi = /(0) in Og . (5.4b) 

(Note that (of )tU=o depends only on uo and not on the values taken by u at times t > 0.) 

Then there exists T G (0, T) depending on uq, f , and Oq, such that there exists a solution 
(v,q) G Z T of the problem fOp . Furthermore, n G C°([0, T]; # 3 (ft£; M 3 ) n ff 3 (Og;R 3 ) n 
ff x (0;: 



Remark 4. In Theorem \13. 6 A assumptions ensuring uniqueness of the solutions are also 
given. 
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Remark 5. If we had not made the notational simplification of Section^ we would have 
to modify \5.ty) by 

[Def Uq N] tSLn = on To, w\ = on dQ, vAuq — Vqo = on Tq, 

[(u[T>efw( Nf + „[ 4,1 (ofof) t (0) + u f ,[ (a*aj) t (0) ] Nj)UU 

= [(c^ kl u f jN j )l 1 U 1 on T , 

and \5. 3b]) would be replaced by q§ = v [Def Uq N] ■ N on Tq. 

Remark 6. The regularity of our solution v G Wt implies that for each t G (0, T] the solid 
domain £l s (t) is of class H 3 . Also, although we have stated our results for three-dimensional 
motion, all of our results hold in the two-dimensional case as well. 

Remark 7. We have stated our results using the convention of Section wherein the 
forcing function f is taken to be defined over the entire domain 0,. It is certainly possible 
to define separate forcing functions for the solid and fluid phase, in which case we would 
need the following regularity: 

f f G L 2 (0,f ;H 1 (n;R 3 )), f ft G L 2 (0, T; L 2 (tl; M 3 )), f fu G L 2 (0, T; HQ n (Q; K 3 )'), 

f s G L 2 (0,f;H 2 (n s ;M. 3 )), f. t G L 2 (0, f; H\n s ; M 3 )), f stt G L 2 (0, f; L 2 {%; R 3 )), 

//(0)65 3 (^;M 3 ) • 

The compatibility condition uAuq — Vqo = on To in h5.l3\) would be replaced by vAuq — 
Vgo + //(0) = / s (0) on T . In the definition of qo, //(0) replaces /(0) and in the definition 
ofw\, /(0) is replaced by //(0) and / s (0) respectively in S7q and CIq. 

Remark 8. Note that the supplementary regularity condition for u(0) and /(0) is due to 
the hyperbolic scaling of the velocity and forcing in the fluid. A parabolic scaling in the 
fluid, which may appear to be more appropriate, would not, however, lead to the necessary 
estimates, except for the case in which the initial solid-fluid interface is flat (which is not 
the case considered herein). This is due to an elastic energy integral (which we shall shortly 
identify) that requires the hyperbolic scaling in order to be estimated. 

Remark 9. Note also the presence of two compatibility conditions for the stresses on Tq, 
which is also a consequence of the hyperbolic scaling. A fluid-fluid interface problem would 
require only one compatibility condition. 

Remark 10. Note that the proof of existence of solutions requires only the "minimal" regu- 
larity assumptions 15. 1\) on the forcing function f ; this is due to our method of proof which 
employs (just as in jS] for the Navier-Stokes equations with surface tension) the Tychonoff 
fixed-point theorem instead of a Banach-type contraction mapping. Note also that unlike 
the case of a free-surface fluid problem, a Banach contraction method does not work for the 
problem that we study herein. We will see later that some additional Lipschitz assumptions 
MS. 3\) are necessary for uniqueness. 

Remark 11. We also remark that our technique is restricted to the case where the elastic 
constitutive law in the solid is either linear or semi-linear. Whereas the paper is written with 
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a linear elasticity law, we can handle in the same fashion and with the exact same methods, 
the case where an extra contribution of the type F(Vn, 77) is added, with F satisfying the 
usual regularity and growth assumptions. In that case the linear problem {7.1]) , defined 
hereafter, which is used in the fixed-point approach would be replaced by a similar problem, 
with {7.1c]) replaced by 

ft r-t ft 



w\ -[<?*' +F{V w, w) = fin (0,T)xU s , 

Jo Jo Jo 

which does not create any additional difficulties with respect to the analysis of the linear 

case. 

The consideration of a quasilinear elastic law such as the nonlinear Saint- Venant Kirch- 
hoff material, involves a smoother functional framework and will be developed in a future 
article. 

6. A BOUNDED CONVEX CLOSED SET OF Wt 

Definition 6.1. Let M > be given. We let Ct(M) denote the subset of Wt consisting 
of elements u £ Wt such that 

\W\\w T < M, (6.1) 

and such that 

with w\ defined in Theorem \5.1l and where we continue to assume that the conditions stated 
in Theorem \5.1\ for the forcing function f and the initial data uq are satisfied. 

Lemma 6.1. There exists Mq > such that Ct(M) is non-empty for M > Mq. Further- 
more, Ct(M) is a convex, bounded and closed subset of XT- 
Proof. We note that if v(t) = u Q + tw\, then v £ C M (T) for M > Mq = \\v\\yy T - The fact 
that Ct(M) is closed follows from Mazur's lemma. □ 

Remark 12. Note also that ifO<T'<T, then Ct'(M) is non-empty. Henceforth, M is 
assumed to be larger that Mq. 

In the remainder of the paper, we will assume that 

< T < f 

where the forcing / is defined on the time interval [0, T]; we will have to choose T sufficiently 
small to ensure existence of solutions to our problem. 

We will need the following series of simple lemmas on the set Ct(M). 

LEMMA 6.2. There exists T £ (0, f) such that for all T £ (0, T ) and for all v £ C T {M), 

the matrix a is well defined, 

and satisfies the estimate (which is independent of v £ Ct(M)) 

ll a llx,ooen T- u2fnf-v>9\\ "I" ll a i 



(0,T;H 2 (n^;K 9 )) + I ' a * I ' (0,T; 1 (S7^' ;R 9 )) + " a " Hl°°(0,T;L 2 (q£ 



+ IKIlL2( 0ir;// 2(^/. R 9)) + ll a ttlli2( 0j T;Hl(n^a 9 )) + H ai *'llL 2 (0,T;L2(Q/ ; ]R9)) ^ C(M) . (6.3) 
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Proof. Notice that in the separable Hilbert space H 3 (Qq ; M 3 ) (for which the Bochner integral 
is well defined), 



77(t) = Id + f 
Jo 



v(s)ds ; 



this together with the Jensen and Cauchy-Schwarz inequalities shows that 

\\V — Id|li<»(0,T;.ff 3 (n£;R 3 )) - ^ ^ I? ,T ' 

and thus 

ll V7 ? - ^L^T-H^m <cVf \\v\\ Xt <cVTVm. (6.4) 

Next, choose R > be such that for any 3x3 matrix b satisfying ||6 — I|| K 9 < R, we have 
detft > |. 

We then see from (|6.4p and the Sobolev inequalities that for T < T = ^rjh, Vr/(t) is 

invertible for t £ [0, T] in Q for any v S Ct(M). From now on, T is assumed to be in 
(0, To). Since 

a(t) = CofV#) in n{ , 

we then see from 1)6. 4JI that 

H a llL-(0,T;^(^;R9)) ^ C '( 1 + V 7 ™) 5 ■ 

Similarly, 

||« - u oll£cc ( o,r; fl 3(n^;]ii3)) ^ C >/TM> ( 6 - 5 ) 

providing 

ll a *lli»(0,T;H 1 (n^R 9 )) - ^ + H U °lli? 3 (^;R 3 ) + ■ 
In the same fashion, 

IK - wi IU«. (0 , W (n^R8)) <CVTM, (6.6) 

providing 

The L 2 in time estimates are established in a straightforward manner from the definition 
of Ct{M), which concludes the proof of the lemma. □ 

Remark 13. Note that Tq also depends on M. 

In the following, T is taken in (0, To) (and M is still taken in (0, Mo)). By the same 
arguments as above, we can easily prove the following results: 

Lemma 6.3. For all v £ Cm{T), we have 

II" " a (°)ll 2 Loc ( o iT;H2( o/ ; R 9)) + K " ^(0)11^(0,^(0/^)) < C ( M ) T ■ ( 6 - 7 ) 

Lemma 6.4. There exists T\ G (0,Tq) which depends on M, and a constant C > which 
depends on uq but does not depend on M , such that we have for all v E Cm{T), 

^^(O^tfSCn^R 3 )) + H U IL°°(0,T;ii" 2 (n/;R 3 )) + H W * ^^(O.T^^oi* 3 )) ~ ^ ' 
The next result concerns potential solid-solid or solid-container collisions for a short time. 
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Lemma 6.5. Let d > denote the infimum of the distances between two distinct connected 
components of (if we have more than one solid in the problem) and of the distance 
between S7g and d£l. Then, there exists T2 £ (0, To) such that for all v £ Cm{T%), 

T2 d 

f T f T 1 

Proof The inequality J \\ v \\ L °°(n f ;R3) - C ^ \ j W v Wh 2 (^m?^ proves the result. □ 

Henceforth, we shall require 

Te(0,T M ), T M = min(T 1 ,T 2 ). 

The series of estimates in SectionslHlandElwill show that M must first be chosen sufficiently 
large, and then T must be chosen sufficiently small. 

The next result is crucial for the derivation of appropriate estimates; while it appears 
that we should require an estimate of qu in L 2 (0, T; L 2 (£Iq ; R)), we are not able to obtain 
such an estimate, and effectively replace it with an estimate of q t in 

Lemma 6.6. For all v £ Cm{T), 

Proof. Let ip(t) = / |ati(i)| 3 + 1 > 1. We then have in the distributional sense 

ij/(t) = 3 f \a tt (t)\ 2 a tu (t) . 



Thus, 

i>'(t) < C ||Ott(t)||^ (n / ;R9) ||Ottt(t)|| L 2 (n / ;R 9 ) > 

which by interpolation yields 

tf(t) < C l|att(i)l| L 3 (n /. M 9 ) ||att(*)|| if i (n /. R 9 ) ||attt(*)ll i 2 ( ^; K 9) 



i.e 



i//(t) < C [^(t)]s l|ot t (t)|| jffl(n /. R9) ||a t t t (i)|| i2(n /. R9) • 
Thus, since ip(t) > 1, 



i>(t) < [if>(0)* + C J q \Mt)\\ Hl{n f.^ ) \\a t tt(t)\\ L2{n f. 
which by ()6.3|) provides 

i/)(t) < [V(0)f +C(M)]I , 
which establishes (j6.1U|) . □ 

Remark 14. Note that in the above L? estimate, the exponent 3 is the limiting case for this 
lemma. 
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Remark 15. Had we not made the notational (constitutive) simplification of Section^ 
we would require the following Korn-type lemma in the Lagrangian setting (this is the only 
mathematical issue that the actual constitutive law hl.le\) requires): 

Lemma 6.7. There exists T% G (0, T) such that for any T G (0, T3) and v G Ct(M), for all 
<f> G H^(n f ;R 3 ) and t G [0,T] 

^(4(t)^+af(t)^J(4(t)^+a 4 fc (t)^) > C IHI^ (n /. R3) • 

Proof. To prove this result, we let a(t) = I + [a(i) — a(0)] and apply (|6.7j) followed by the 
Korn inequality. □ 



7. The basic linear problem 



Suppose that M > Mo, T G (0,Tm) and u G Cr(M) are given. Let i] = Id+ / i> and 

Jo 

let be the quantity associated with r\ through ()3.1I) . 

We are concerned with the following time-dependent linear problem, whose fixed-point 
w = v provides a solution to ()3.2|) : 

w\ - v{a\a\w\ k ),j +{a k q) lk = f o rj 



a k w\ k = 



Jjkl I „„k 



(0,T) x^, 
(0,T) xfi£, 


(7.1a) 
(7.1b) 


(o,T)xng, 


(7.1c) 


(o,r) xr , 


(7.1d) 


in (0,T), 
ft x {t = 0} , 

n x {t = 0} , 


(7.1e) 
(7.1f) 
(7-lg) 





vw\u afajNj - qa{Nj 

w(t,-) G H^n-m 3 ) 

W = Uq 

rj = Id 

The following regularity result will be of paramount importance in our analysis: 

Theorem 7.1. Given f and uq satisfying the assumptions of Theorem I5.il there exists 
M > 0, T > 0, such that for any v G Ct{M), there exists a unique solution (w,p) G Zj- of 
f71\ >. Furthermore, w G C T (M). 

Sections and EH are devoted to the proof of this theorem. In the following, we set 

N(u , f) 2 =(1 + IKI|^ S(n /. R3) + Kllfp^RS) + ll/(0)llfr3 ( n;R3) + ll/llia( ,T;iF»(n 
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8. Preliminary results 

8.1. Divergence, extension and regularization type results. We first state the fol- 
lowing result, whose proof follows the same argument as for the case of a smooth boundary, 
with the exception that the regularity results for elliptic systems of ^2] are used instead of 
the more classical results wherein the boundary is smooth. 

Lemma 8.1. Let O' be a domain of class H k (k > 3). 

Then, for < m < k — 2, there exists a continuous linear operator 

L{n'):{{d,r)GH m (Q';R)xH m+0 - 5 {dn';M: i )\ [d=[ r ■ n} -> H m+1 (Q'; R 3 ) 

Jn Jan 

such that u = L(Q')(d,r) satisfies 

divu = d in £1 , 
u = r on dfl' . 

Furthermore, the operator norm of L(Q!) remains bounded as the norm of the charts defining 
Q! stays in a bounded set of H k . 

We will need the following extension 

Lemma 8.2. Recalling that SIq is of class H 3 , for each 1 < m < 3, there exists a continuous 
linear operator 

E : H m (n f ;R 3 ) D fl^ n (n£;R 3 ) -> H m (fl; R 3 ) D H%(Q; M 3 ) 
such that E(u) = u in Oq. . 

Proof. The result is well-known in the case that S7q = W\_; let n : R 3 + -> K 3 denote this 
extension operator, and let {^i}^ denote a collection of charts in a neighborhood of Tq 
(each is a map of class H 3 from the unit ball in M 3 into an open set containing a 
coordinate patch of To), and let {6i)f =l denote the associated partition of the unity. We see 
that 

N 

F{u) =^II[(0 i u)o^]o#r 1 
i=i 

is an extension of u into a neighborhood of Tq. By introducing a smooth cut-off function £, 
equal to 1 in and equal to in the complementary part of this neighborhood included 
in Oq, we see that E(u) = £ F(u) satisfies the statement of the lemma. □ 

In a similar fashion, we can also extend from Q to IR 3 , with the same arguments. 

Lemma 8.3. There exists a linear and continuous operator E g from H m (Q; M 3 ) into i/ m (R 3 ; M 3 ) 
(for each 1 < m < 3) such that E g (u) = u in Q. 

We also need a regularization lemma for the coefficients a\ , which we shall use to obtain 
estimates for the solutions of the regularized problems (whose coefficients by definition use 
these regularized coefficients) ; we will then pass to the limit as the regularization parameter 
tends to zero. 
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Lemma 8.4. Let v G Ct(M) and rj = Id + J v. Then, there exists a sequence v n G 

Vp 9 (T) = {ue L 2 (0,T;H 3 (n f ;R 3 ))\ u t G L 2 (0, T; H 3 (Q f ; R 3 )), u tt G L 2 (0,T; H 3 (Q f -R 3 ))}, 
such that v n (0)\ f = u \ a f, v nt (0)\ of =wi\ f, and 

K - «||y*(T) + ||K - W )«IIl-(0,T;L2(^;M3)) ~» 0. 

Proof. Let p G V(B(0, 1)) be such that p > and / p = 1, and let p n (x) = n 3 p(xn) 

ifl(0,l) 

denote the usual mollifier. 

From Lemma IO for any t G [0,T], let v(t) = E(v(t)), so that v G v/(T) n Kf(T) with 
ll^llv^T) + ll^lly s 3 (T) < C IMIv^t) • We extend to R 3 by setting v' = E g (v). 

Then let v n be defined for any t G [0, T] by 

v n (t) = p n *v r (t). 

From the properties of the space convolution, we have that v n G Vj e9 (T) n V 3 (T) and that 

11% - V'\\ V 3(T) + \\Vn ~ ^'Wv^T) + II (v n ~ v') tt || L oo ( ,T;L 2 (c£;R3)) ~* aS n 00 • Tllis in tum 

implies that ||u n — i>||y3(y) + ||(f n — ll^oo^ T-L 2 (n / -R 3 )) ~^ ^ as n ~~ * 00 • Now, for the initial 
conditions, let us define 

v n (t) =u + twi- / (t' - t){v n )tt dt! . 
Jo 

(the Bochner integral being well-defined in the Hilbert space H 3 (Q f ;R 3 )). We then have 
that v n (0) = uq, v nt (0) = w\. 
Moreover, 

v n = v n + E g (E(u )) -p n * E g (E{u )) + t [E g {E(w x )) - Pn * E g (E(wi))), 
which yields \\v n - v\\ V 3^ + \\(v n - v)tt\\ L oo( QiT . L 2( n f. R 3)) -> 0, as n -> 00. □ 
Remark 16. Our construction does not necessarily yield v n = on (90. Consequently, with 

r\ n = Id+ v n we do not have that r/ n (0) = Q. It, nevertheless, does not matter for the 
Jo 

purpose of our analysis. 

Remark 17. In the following, we will solve {7.1}) as the limit as n — > 00 of the solu- 
tions w n to the problems {7.1]/ associated to these regularized v n . The interest of this 
regularizing process is that for a given n, a(j] n ) and its first and second time derivatives 
are in L°°(0, T; H 2 (Qq;R 9 )) C L°°(0, T; L°°(f^; R 9 )) and its third time derivative is in 
L 2 (0,T;H 2 (n f ;R 9 )) C L 2 (0, T; L°°(ftjJ; R 9 )) which is necessary in order to get the exis- 
tence of regular solutions to {7.1] ). These bounds in those spaces of course blow up as 
n — > 00 (except for the estimate for a(rj n ) ). 

Nevertheless, using the fact that \\v n — v\\ V 3^ — > as n — > 00, a(rj n ), and its first, 

second and third time derivatives satisfy the same type of estimates as \b\ \6.1U\) and 
{6. 7\) , respectively, with a constant C{M) which does not depend on n. This fact will be 
used, together with interpolation inequalities (that hold since the solutions w n are regular) 
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in order to get estimates in Yt for w n which are independent of n. By weak convergence, 
this will provide our smooth solution to \7.1\) . 

We will also use the convention of denoting the regularized velocity fields v n by v, and 
the corresponding regularized matrix a{ii n ) by a. 

Remark 18. Since the fluid forcing in \7.1\ ) is given for], we need to extend f to M 3 . 
Hence, when we solve this problem with the regularized coefficients arising from v, we in 
fact implicitly use the extension E g (f). This extension has the same regularity as f with 
M 3 replacing £1; this follows from the fact that E g commutes with the time derivative. For 
notational convenience, we shall continue to denote the extended forcing function by the 
same letter f . 

8.2. Pressure as a Lagrange multiplier. 

LEMMA 8.5. For all p G L 2 (Q,q;R), t G [0,T], there exists a constant C > and 4> G 
-fT^rijM 3 ) such that aj(t)(f> l ,j = p in and 

<C||p||2 (8.1) 



lifi(fi;R 3 ) ^ ^llL2 (n /;i 

Proof. Let p\ G L 2 (Oq; R) be such that p\ > in f2g. Let p be defined by p = p in f2jj and 
J Qf pdeWr,(t) 

Pi detV7?(t) 

po n(t)~ l dx = p(x) detV rj(t)dx = , 
T)(t,n) Jn 

we then see that <t> = L(rj(t, fi))(p o r/(i) _1 , 0) o 77(f) G i^(ft;M 3 ) satisfies 

div(</> o 77(f) 1 ) = p o n(t)~ in n(t, Q) , 

and thus 

ai(t)4> % ,j = dw(cf) o ^(i) -1 ) o 77(f) = p = p in fijj . 

The inequality (|8.1() is then a simple consequence of the properties of L and of the 
condition v G C T {M). □ 

We can now follow |21j . Define the linear functional on Hq(£1; R 3 ) by (p, a^(t)tp^,j )^2^q/.^j 
where tp G Hq (Q; R 3 ). By the Riesz representation theorem, there is a bounded linear op- 
erator Q{t) : L 2 (n f ;R 3 ) -> H^(Q;R 3 ) such that 

\/tp G Hq(£1; R 3 ), (p, cti(*)^>i)L2 (n /. R) = (Q(t)p, <p) H ^(U^)- 
Letting tp = Q(t)p shows that 

\\Q(t)p\\Hu^m < c \\p\\mnfm ( 8 - 2 ) 

for some constant C > 0. Using Lemma 18.51 we have the estimate 

I^"l2(^ ; r) - c WQ( t )p\\H^nm^H^(nm < c 'll < 5( i >lli? i(n;R3)lblU2 (n /. R) , (8.3) 
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which shows that R(Q(t)) is closed in H^(n-,R 3 ). Since V v (t) C i?(Q(t)) ± and ^(Q^)) 1 - C 
V v (t), it follows that 

H^R 3 ) = R(Q(t)) ®Hi (am V v (t). (8.4) 
We can now introduce our Lagrange multiplier 

Lemma 8.6. Let £(t) G i2" _1 (fi;M 3 ) 6e such that £(t)cp = for any ip G V«(t). Then there 
exists a unique q(t) G L 2 (f2Q;R), which is termed the pressure function, satisfying 

VveH^M 3 ), 2(t)(<p) = (q(t), ay\ 3 ) LHn f. m - 

Moreover, there is a C > (which does not depend on t £ [0,T] and on the choice of 
v G Cm(T)) such that 

Proof. By the decomposition (|8.4j) . for u G ifo(f2,M 3 ), we let ip = V\ + t>2, where t>i G V„(i) 
and V2 G R(Q(t)). Prom our assumption, it follows that 

£(t)((p) = £(t)(v 2 ) = (^(i),«2)ifi(n,R3) = W>(t)> ( P)H%(njP)> 

for a unique ip(t) G R(Q(t)). 

From the definition of Q(t) we then get the existence of a unique q(t) G L 2 (Qq;R) such 
that 

V^Gtf^lR 3 ), 2(t)(cp) = (q(t), 4^)^ {a f;Ry 
The estimate stated in the lemma is then a simple consequence of ()8.3[) . □ 

We will also need a version of the Lagrange multiplier lemma for the case where £,(t) G 
1 (SIq ; R 3 ), which implies an estimate on the pressure modulo a constant. We first have 

LEMMA 8.7. For all p G L 2 (Ql;R) such that / p det Vn = 0, t G [0,T], i/iere exists a 
constant C > and ^ G i^Q (J1q ; M 3 ) suc/i i/ia£ a^t)^ 1 ,j = p in $7q and 



Proof. Since 

/„(/ o,;,! 



/ 



po r/(t) 1 dx = I p(x) detVn(t)dx = 



we then define 4> = L(r/(t, OjJ))(p o 77(f) ^0) o ^(t) G i?o(0^; 

The inequality (|8.1I) is then a simple consequence of the properties of L and of the 
condition v G C T (M). □ 

In a completely similar fashion to the proof of Lemma 18.61 we can now establish our 
second Lagrange multiplier 
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Lemma 8.8. Let £(t) G ff -1 (^;M 3 ) be such that £(t)tp = for any tp G V„(t)nH(J(fi£; K 3 ). 

Tfterc i« e » s(s a u« iq ue „( t ) 6 £"(ng;R), .Oft* / V, = 0, W kich termed 

the pressure function, satisfying 

Vc^G^^M 3 ), £(t)( v ) = (<jr(t), afc^ ) i2(n /. M) - 

Moreover, there is a C > (which does not depend on t G [0,T] and on f/ie choice of 
v G Cm(T)) such that 

Remark 19. Tae /ottr previous lemmas do not rely on the fact that v = on d£l. Therefore, 
they are also true for the case where the coefficients a are associated to v. The important 
point is that the estimates i|6'.5)) . \b\ 1U\) and j6'.7| ) are also satisfied by the regularized matrix 
a and velocity v. 



9. Estimates for (|7.1j) : the case of the regularized coefficients 

9.1. Weak solutions. 

Definition 9.1. A vector w G V v ([0,T]) mfre G V v (t)' for a.e. t G (0, T) is a u>eafc 
solution of i7.1p provided that a.e. t G (0, T) ; 

(i) (tot, $ + ^Kw>V , al<P\s ) L 2 (n f. R 9 ) + (c ijM f w k d , 4>\j )i2 (n g;R) (9.1a) 

«* 

= (^'^^(n^RS) + (/,0)i 2 (ng;]R3), G V v (i) , and 

(ii) w(0, •) = «o, (9.1b) 
/or a.e. < i < T , where (•, •) denotes the duality product between V v (t) and its dual V v (t)' . 

9.2. Penalized problems. Whereas the existence of a weak solution can be proved directly 
in the space V„([0, T]), with wt G V v (t)' , this framework is not amenable to finding the 
pressure estimate required by our analysis. (Even for the well-studied problem of the Navier- 
Stokes equations on a fixed and smooth bounded domain, the weak solution only provides a 
pressure estimate of the form f Q p G L 2 (0, T; L 2 (Qq; K)).) A penalized form of the problem, 
however, together with the penalized form for the time-differentiated problem, provide the 
correct pressure estimate in the limit as the penalization parameter tends to zero. 

As we noted following Lemma 18.41 we will work with a regularized sequence of velocities 
v n , and we shall generically denote elements of this sequence simply as v, and the associated 
regularized matrices a(%) as a. 

Given the regularity assumptions in (fS"Tj) . f t G C([0, T];L 2 (£l; M 3 )), so that f t (0) G 
L 2 (n f ;R s )). 

Then, let w 2 G L 2 (0;M 3 ) be defined by 

w\ = uAw\ + u((ajaf) t (0)ul k ), j +F t (0) - {(4) t (0)q ) :j -q 1:i in n f , (9.2a) 
wi = fm + [c m u k ,l],j in fig, (9.2b) 
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where q\ 6 H 1 (f^jR) is defined by 

A Q1 = ^AKK + (F%(0) + K(a>f) t (0H, fc )„ -((^) t (0)g )„ ] 

+ 2(^(0)™*,,- +(aj) tt (0)4, j in fi£ , (9.3a) 
= u[V Wl N-N + (afaj) t (0)ul k NjNi] - c ijhl 4,i NjNi on T 

+ 90 K)t(0)iV J ^ (9.3b) 

|| =F t (0) • iV - [(af )t(0)g ],j +isA Wl ■ N + is((aj of ) t (0)u , jt )„■ iV, on 90 . (9.3c) 

Once again, we remind the reader that (a^)t(O) and (a^)it(O) depend only on no and w±, 
and we note that they are equal to (a])t(0) and (a])tt(0), respectively. 

Letting e > denote the penalization parameter, we define w e G W([0,T]) to be the 
unique weak solution of the problem (whose existence can be obtained via a standard 
Galerkin method in a basis of .^(fi; R 3 )): 

r-t 



(i) (w et , (f>) + v(a r k w l e , r , a s k <p\ s ) £ a (n /. R) + (c lJ J w € , ^(ngjR) 

+ (^n?^ -g - %, «U fe >« ) L 2 (Q / ;]R) = (i 7 , 0) L 2 (n / ;M3) + (/, ^)i2(Qg;R3), (9.4a) 
V</> G i^^R 3 ), and 

(ii) n>(0, •) = n , (9.4b) 

where (■, •) denotes the duality product between ^(OjR 3 ) and its dual. 

9.3. Weak solutions for the penalized problem. The aim of this section is to establish 
the existence of w e , as well as the energy equalities satisfied by w e and w et , and the energy 
inequality satisfied by w ett . It turns out that the exposition is simplified if we first study 
the twice differentiated-in-time problem, that we introduce now. 
Step 1. Galerkin sequence. 

By introducing a basis (ei)fZ\ of Hq(Q; R 3 ) and L 2 (VL; R 3 ), and taking the approximation 

i 

at rank I > 2 under the form wi(t, x) = ^ Vk{t) e k (x) , and satisfying on [0, T] 

k=i 

(i) (wittt, 0)l2(^ ; r3) + v(a\wi tUr , a s k (j), s ) £2(n /. R 3) + (c ljkl wit,i , 4>\j ) L 2 (n g. M) 

+ v((a k a k ) tt wi,r , <P,s ) L 2 (n /. M3) + 1v{{a r k a s k ) t wi Ur , & a ) £2(n /. R3) 

- ({afqih, 4>\j) L 2 {n f. R) 

= ( F tt,<P) L 2( n f. R 3) + (/tt^)L2(Qg;R3), V0 G span(e!,...,ei) , 

(ii) n; U ((0) = (n; 2 );, n> u (0) = (n>i);, n>;(0) = (n )z, in O , 

where qi = qo + t q\ &i w hj an d (^2)/ denotes the L 2 (S7;R 3 ) projection of n>2 onto 

span(ei, ei), and (n^i); and (no); denote the respective ^(OjR 3 ) projections of w\ and 
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uq on span(e\, ...,ei), we see that the Cauchy-Lipschitz theorem gives us the local well- 
posedness for wi. The use of the test function (wi)u in this system of ODEs (which is 
allowed as it belongs to span(e.\, ei)) gives us in turn the energy law 

^ ^ II I|2 ( ~T ~s \ 

2 ~dt * llz, 2 (^;]R3) + v \ a k w lttir , a k W ltt>s) L 2(Qf. R 3) 

+ ^(c ljkl u>h ,i ,wti,j)Li(n°-,R) + v((a r k a s k )ttWi,r , witt,s) L 2 {n f. R3) 
+ 2is((a r k a s k ) t w lt , r , w ltt , s ) L 2 (n /. R s) - {{a\qi)tt, w Ut,j ) L 2 (n f. R) 

After transforming the term with {qi)tt (since it involves Vwi tt ) and integrating this 
relation from to t G (0,T), 



1 f 1 

2 INtt Iliads) + ^ y (al w ltur, a s k w ltt , s ) L2 ( n f. R3) + -(c* jM ^i , )L2(n-K) 

+ e / f f Qltt+ f I . Qltt P(rt)tmij +(al) tt wi\ j }-2 f f {a^t^w^j 
Jo Jn(, Jo J% Jo J% 

-u 

) ii; Ut) L 2 ( -^/. R 3) + y (ftt,Wl tt ) L 2( n s. R ;iy (9.5) 

By noticing that, e being fixed, the fourth term of the left-hand side of this inequality 
involving the square of (qi)tt acts as a viscous energy term, and taking into account the 
L oo (0,T;L oo (^;M)) bound of each one of the regularized coefficients a\ and their first and 
second time derivatives, we then get 

l\\witt(t)\\l2 (Qm + V -=\ \\VwiuWl f + hj^wtli (t),wi\,j (t)) L 2 i%m 



2 II ™vv \\L*(\l;K?) ' 2 J Q 11 ' wlIIII L2(^;IR9) ' 2 

rt r rt ft' ,■'< rt' 



< C N(u ,f) 2 + (Ftt,w ltt ) L2{n f. R3) + J (ftt, w Ht ) L 2 {n s. R3) , 

where C e depends on the regularizing parameter of a and on e, but not on I. By GronwalPs 
inequality, we then get an estimate on each of the integral terms multiplying C e which in 
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turn implies 

\\ w itt(t)\\ 2 L 2 {nm + J ll Vw mllJ 2(n / ;R9) + 2^ w it* V)>i»iU (t))LHn°;R) 

+ 6 / [ qi 2 tt <C e N(u J) 2 . 

Jo Jni 



Step 2. Weak solution w e of the penalized problem and its time differentiated 
problem. 

We can then infer that wi is defined on [0, T], and that there is a subsequence, still 
denoted with the subscript I, satisfying 

Wi^We inL 2 (0,T;^(O;M 3 )) (9.6a) 

w lt ^w et mL 2 (0,T;H^(n;R 3 )) (9.6b) 

witt^wett in L 2 (0,T;L 2 (ft;]R 3 )) and in L 2 (0, T; F 1 ^; M 3 )) (9.6c) 

g /tt - g ett in L 2 (0,T;L 2 (^;M 3 )) , (9.6d) 

where 

q e = q + % - -a{we,j ■ (9.7) 

From the standard procedure for weak solutions, we can now infer from these weak con- 
vergences and the definition of wi that w ettt G L 2 (0, T; H~ 1 (Q; M 3 )). In turn, u> e44 E 
COQO^IjF-^^R 3 )),^ G C7°([0,T];L 2 (^;M 3 )),u; e G C°([0, T]; M 3 )), with w e (0) = 

uo, w €t (0) = wi, w ett (0) = w 2 . 
We moreover have for w\ t 

(i) (wut, 0)l2(h ; m3) + v{a r k wi t , r , a|(/>, s ) L 2 (n /. K 3) + {c l3kl wi k ,i , 0*y )l,2(ng;M) 

+ v((a r k a s k ) t w hr , ^ s ) L 2 (n /. K3) - ((aUih 4>\j) L 2 {Q f. R) 

= (Ft, 0) L 2 (n /. R 3) + (/t) 0)l 2 (^ ; ir3) + q(0), M4> E span(ei, e/) , 

(ii) wi t (0) = (wi)i, wi(0) = (uo)i, in O , 

where ci(4>) G M is given by 
ci(4>) = ((w 2 )/,0)l2(c ; k3) + v(a r k(°) w u,r , ai{0)(l), s ) L2{n f. R3) + (c tjkl (wi) k ,i,(l)\ j ) L 2 {n s. R) 
+ v((flla s k ) t (Q)(wi)i„ , ^ s ) Z 2 (n /. R 3) - ((alqi)t(O), 4>\j )^(n f -M.) 

Thus, q(</>) converges to the same expression, where the approximate initial data (wi)i are 
replaced by the actual initial data Wi (i = 0, 1, 2). From our compatibility conditions Q5.2JI 
together with 1)9.2(1 . this leads us to 

l|Q||/r-i(n;K 3 ) °> as / ^ oo . (9.8) 
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Similarly, for wi 

ijkl I „,,k ii 



W £2(fi. R 3) +v(a V k Wl,r , a l^>s)£2(fi/;]R3) + (C* J J Wl ,1 ,<f>\j ) L 2(nS;M.) 

= (F, ^(q/.rs) + (/> 0)i2 ( ng;R3) + ci(4>)t + G span(e 1: e z ) , (9.9a) 

(ii) u>,(0) = (u )i, in U , (9.9b) 
where d/(^>) G M is given by 

<W) = ((wi)/,0)l2(q. R 3) + 1/(0^(0)1X0^ r , 5 K O )0>s) L 2 (n /. R 3) 

Similarly as for ci(4>), from our compatibility conditions Q5.2|) 

INIh-i(c ; r3) -> 0, as / ^ 00 . (9.10) 
We can thus infer now that at the limit u> e satisfies for all G L 2 (0,T; Hq(Q;R 3 )), 



Wf.t,4>)l?{n-$i?) dt + V (d k W e , r , o|0, s ) i 2( n /. 



(-F, <A)l2(^;IR3) + (/) 0)L 2 (Cg;R3) (ft , (9.11) 

which, combined with io e (0) = uq, shows us that w t is a weak solution of (|9.4[) . 
Moreover, u; ei satisfies for all G L 2 (0, T; //q (J7; ] 



(^ett,<^)x2(n ; K3) dt + u J ((a s k a r k w e , r )t, 0>s)i2( n / ;K 3) * 

/•T /-T 

+ / (c^ fc 'u; e fc , i ,^, i )L2(cg;R) eft- / ((4fc)t> <l> d) LHu f. M 



T 

(F t , 0) i2(n /. R 3) + (/t,0)L2 ( ng ;K 3) . (9.12) 

Step 3. Strong convergence for the Galerkin approximation. 

Since w e G L 2 (0, T; //q (S7; M 3 )), we can use it as a test function in l|9.11|) . which provides 
us on (0, T) with the equality 
1 rt 

2ll W e(*)lli 2 (n;K 3 ) + U J Q (flk w e,r, d s k W e , s ) L 2 (n / ;R 3) dt 

1 /"* 

= ^H^lli 2 ^ 3 ) + J (- F '^)l2(^ ; ]R3) + (f,w e ) L 2 {n s. R3) dt . (9.13) 



2-1 
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Similarly since wi(t) E span{e\, ...,ei) for all t G [0, T], we can use it as a test function in 
(|9.9|) . which gives us 



1 2 rt 

= 2H( U °)^i 2 (0;R3) + / ( F > w *) + (/'^)L2(Cg;R3) dt 
J 

+ / tdi(wi) +ct(wi) dt. (9.14) 

JO 

By (HSU) , (|9~H|) and (|9~TU|) . we then infer by comparing (j9~H|) and lfD~T3|). that as I -> oo, for 

all t G [0,T], 

1 



2 H^(*)llz,2(Q;R3) +U / (aJ>J, r , d S k Wi, s ) L2{n f. R3) dt 
— J 

1 /"* 



2 



1 



o 

/ ,„ k , / * . //+ / \\n II 2 



+ 2 (c 7o We ''' 7o V ' i)i2 ^« * + e / ll<?e|l ^«^ ' 

which gives in turn the strong convergences 

w t ^w e inL^TjJET^njkR 3 )) (9.15a) 

wi^w e in L 2 (0,T;L 2 (ft;IR 3 )) (9.15b) 

f wi^ f w e in ^(O.rjff 1 ^;! 3 )) (9.15c) 
Jo JO 

gj g e in L 2 (0, T; L 2 (^; R)) . (9.15d) 

Since u) et G £ 2 (0, T; Hq (fl; R 3 )), we can prove in a similar fashion the strong convergences 

wi t -» ^et in i 2 (0, T; H 1 ^; R 3 )) (9.16a) 

w /t ^w £t in L 2 (0,T;L 2 (O;M 3 )) (9.16b) 

u; z -► w e in L^T^^R 3 )) (9.16c) 

g It - g £i in L 2 (0, T; L 2 (5^; M)) . (9.16d) 

Step 4. Energy inequality for w tt f 
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By using the relation 

J- ,, ,,9 / /_». _ e , 1 



+ e / f t Qltt+ f [ . Qiu [ 2 («i)^t>j ] " 2 / [ f {ai) t qi t wiit,j 

Jo Jcii Jo Jnl Jo Jnl 



X^ihrnmluj + v j ((a r k a s k )ttwi, r , wzttj«) £ 2( n /.RS) 

'0 ^° 

ft 

+ 2z/ / ((o^a|)fu; u , r , w ett , s ) L2(nf , 



o 



'L 2 (QJ ;R3) 



< C N(u , f) 2 + J (F tt 

> W ltt) L 2( n f ;R3) + / (ftt> Witt, 



L 2 (Cg;R 3 ); 



and the weak convergences (|9.6|) . along with the strong convergences ()9.15(l and (|9.16f) . we 
then get 

1 /"* 1 

2W w ttt\\h(n;K3) + v \ {a r k w ett,r , a s k w ett , s ) L 2( n /. R 3) + 2^ c%j We t^ > We t>j )i 2 (^;iR) 



+ e / [ q e tt+ [ f qea [2(4)tW £ lj + (aJ)ttW e i ,j ] - 2 I [ m^w^ 
Jo Jnl Jo Jnl Jo inj 

r-t 



(al)ttqeWel t ,j +p / {(a r k a%)ttw e , r , w ett , s ) 2( f 3) 
'o Jfi f Jo < o> J 

<C N(u ,f) 2 + ^ (F tt 

9.4. Existence of tu, fi^, u5#, uniqueness. In this section, we establish the existence of 
w, and its first and second time derivatives by taking the limit e — ► 0. The inequality 
(|9.53|) proved at the end of this section, holds for any regularized velocity field v = v n , 
independently of n, and requires in its proof, strong convergence results from their penalized 
counterparts since the regularity that we take on the data does not allow us to view ibtt as 
a weak solution of a variational problem. 

Theorem 9.1. Suppose that uq and f satisfy the conditions stated in Theorem \5.1i Then, 
there exists a weak solution w to the problem \7.1]/ with the mollified coefficients replacing 
the actual coefficients. Moreover, w is in L 2 (0, T; Va(-)) and is unique, and vbt £ W([0,T]). 

Proof. Step 1. The limit as e — > 0. 

Let e = — ; we first pass to the weak limit as m — > oo. The energy law (|9.13|) shows that 
there exists a subsequence {wj_} such that 

w^^w in W([0,T]). (9.18) 
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Moreover, since (|9.13|) also shows that ||a^uA_,j || L 2( T-L 2 (Q f -M^ ~^ ® as 771 ~^ °°> we then 



have || ajU) J ,j II L 2( 0jT;L 2(q/ ;R )) — 0, i.e. 



w G Vfi([0,T]) . (9.19) 

Step 2. The penalized time differentiated problems and estimates independent 
of e. 

Thanks to (l9~THl) and (l9~T71) . we have w ef £ L 2 (0, T; M 3 )) n C°([0, T]; L 2 (ft; M 3 )). 

We can thus use it as a test function in ()9.12j) . which gives us for a.e. t G (0, T) 

1 /•* i 

2 11^*11^2(^.^3) + v I (a k w et , r , a k w et , s ) L2{n f. R3) + -{c J w e ,i ,w e ,j J^n^R) 

J 

JO 7f7^, JO Jf2J JO JjJj JO Jf2£ 

< C lkl|||a(n;R3) + J ( F tl W et) L 2 in f. R 3) + (ft,W tt ) L 2 {n s. R:i) . 

At this stage, we remove the time derivative from the q tt term in this inequality by inte- 
grating by parts: 

qet(fi)tU>e,j = ~ I Qe{{4)t w e\j )t + (qe{4)t w e,j )(*) - <R)(a?)t(0Ky 
JO 

(q e (0) = qo by divuo = 0); we then infer by the regularity of a and of w t that 

1 ll n2 ^ f / ~v ~s \ ^ / iihl k i \ 

2\\ W tt\\L2(n;M?) + 2 l n \ a k w et,r, a k W et , s ) L 2^f. R 3) + ^(c W e ,| , W e! j Ji,2(fig; K ) 

+ C N(u J) 2 + C J \\w et \\ 2 L 2 (nm , 



< C J \\qef L2{n f ;R) +6 he(t)f L2{n f ;M) +C C S \\VW £ ^ ' 



where 5 > is arbitrary, and C denotes a generic constant depending on the smoothing 
parameter n implicit in a. 

Note that it is the presence of 5# which requires the use of the regularized coefficient 
matrix a; this is due to the fact that a#(i) is not in L°° as the presence of Vw e and q t 
(both taken in L 2 ) would require. In order for us to be able to obtain consistent estimates 
later on which are independent of the regularization process, we must require the pressure 
of the penalized problem to be in H 1 (for a.e. t G (0,T)); this requires difference quotient 
methods. In order to achieve this, we first define this pressure function to be in I? (a.e. 
t G (0, T)), and then find estimates which are independent of the regularization of a. Thus, 
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in (0,T), 

2 



ll^t( t )llL 2 (Q;R3) + l|w e (i)||tfl(ng ; R3) + J H^Htfl^RS) 
^ 6 l H^H^ (n /. R) +* he(t)\\ 2 L2(n f m 

+ CC 5 \\VwS)\\% {n f m + CN( Uo ,f) 2 + C^ |K t ||| 2(n . R3) • (9.20) 
By the Lagrange multiplier Lemma 18, 61 we also have 



2 

+ N(u ,f) 2 } , 

which coupled with (|9.20|) and (|9. 13f) . gives for a choice of S > small enough 

t 

2 ^ r / ll~ I|2 i Ilr7„„ / j. \ 1 1 2 i at/™. -f\2i 



ft 

+ C I \\Wet\\ L 2^ n . R 3) ■ 



Since / ||Vt« e (t)||^ 2 / „ < C N(uq, f) 2 , we get by GronwalPs inequality an estimate on 

Jo L (il ;K ) 

/ ? f L„ n / B , which in turn provides 

< C[||V^(t)||^ (n , ;R9) +Ar(uo,/) 2 ]+C^V^IIi 2( n ; K3) • (9-21) 



2 



Combined with (|9.2U|) . still for 5 > small enough, this also gives 

< CN( U0 jf + C |ke(t)||^ (n /. R3) + C'^Ve t |l! 2 (0 ;R 3 ) 

By Gronwall and ()9.13|) . we first deduce a bound on y ll^etll^^.^) which in turn provides 
us with 

2 , „... , M „2 , / „... i,2 . CN(U J) 2 

+ C\\w e (t)\\lunf.^- (9-22) 



^i(*)lli2(n;R3) + IK(*)llffi(ng;R3) + ^ IKtH Hl(n /. R3) - C N ( u ^f) 

H 1 ^ 



Step 3. An estimate of w €t on [0, T] which is independent of e. 
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rt 
/0 



By using w e (t) = uq + w et , we see that 

Jo 



\\ w £t(t)\\L2(n-,m.3) + \\ w t\\m(n°-,m.3) + / 11^*11^1(0^; 



<cN(u , f y + Cl t i [l^tll^^ + ciitioii^^ 



<cN(u , f y + Cl t y o iK*||^ (n / ;R3) , 

where we denote by C\ a constant, dependent on the smoothing parameter of o (but not 
on e), which will remain unchanged in the following estimates. 
Now, we see that for any < t < t\ = Min(T/v/ , we have 

l f l 

\\ W et(t)\\ 2 L 2 inm + \\We\\ H l in s. R3) + ~ \\ W tt\\ 2 H1{Q f. R3) < C N(U J) 2 , 



rti 

which with w e (ti) = u + / w et gives 

Jo 



\ w ^\\m ( nfm- d N{uoJ)2 ■ (9 - 23) 



Next, we take t > t\ and let w € (t) = w € (ti) + / w et ; we have from 1)9.22(1 and (|9.23|) that 

Ju 



lk£t(*)lli2 ( n;R3) + |k e || H l (nS ;RS) + J \ \ W et \ \ Rl {Qj f . r3 ) 

< C N(u , ff + c l{ t- h) f IK t ||^ ;R3) + c \\Mh)f HHQ f. m 



<C iV(uo,/r + d (i-tx) / IKtH^/.j 







Now, we see that for any t\ < t < 2 t±, we have 



1 /"* 

lket(*)|lL2 ( n;RS) + IMIffi(^;R3) + / II II ff l (n / ;R 3) < ^K, /) 2 , 



/•2ti 

which with w € (2ti) = uq + w et gives 



11^(^)11^(^3) < C N(u jf. 
We then see by an easy induction argument that for any t E (0, T), 

Kt(*)ll£»(n;R3) + IMlWg;R3) + ~ f Q IMI^RS) < ^K, /)' • (9.24) 
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(As is evident in the proof, the constant C grows as T increases and thus depends on T.) 
Thus with (pHTJl . for all t G [0,T], 

\\«M\l 2(n f m <CN(u J) 2 . (9.25) 

Step 4. Weak convergence and limit problem. 

Since w e also satisfies (|9.13jl . we thus deduce that for the choice e = ^- there is a 

subsequence, still noted w j_ , such that 

m i 

WJL _ -» w in L^O^iT^tyR 3 )) (9.26a) 

w_L in L 2 (0,T;L 2 (ft;R 3 )) (9.26b) 

q_L~^q in L 2 (0,T;L 2 (^;R)) . (9.26c) 

By the weak convergences ()9.26j) . we infer from ()9.11() that at the limit, for each <j) G 

T r-T 

(w t ,(f>)L2(n;R3) dt + v J (a r k w, r , a|0, s ) i2(n / ;B , 3 ) dt 
+ [ (c ijkl * fc ,z,^,j)L2(ng;R) dt- (q, a l k k ,i) L2{n f. R) dt 

( F ' ^) L2(^;R3) + (/' 0)L 2 (f}g;K 3 ) ^ • ( 9 - 27 ) 

Now for the initial condition, we notice that w G C°([0, T]; L 2 (ft; R 3 )). From the following 
identities which hold in L 2 (f2;] 



w(t) = W(0) + / W t , W e (t) =U + Wet, 

Jo Jo 



we deduce from the weak convergence of / w t t to / wt in L 2 (0, T; L 2 (f2; R 3 )) that w(0) = 

J o Jo 

uq m L 2 (n;R 3 ). Combined with J07J), this shows that w is a weak solution of ()9.1|) 
associated to v. 

Now, let us prove that the sequences in (|9.26|) in fact converge strongly. 
Step 5. Strong convergence. 

Since w G L 2 (0, T; H^(Q; R 3 )) we can use w as a test function in ()9.27|) . which provides 
an energy law that we can compare to (|9.b-{|) . By using the weak convergence in (|9.26j) . 
and the fact that \\a l k w e k ,i || i2 ( t-l 2 (c / -r 3 )) ~~ > ® as 6 ~^ ® from (|9.13j) . we deduce from this 
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comparison that for any t G [0,T], as e — > 0, 
1 /"* 







^\\w(t)\\ L 2 {n . M3) + iy / (a^,r, a|w, s ) L2(Q /. R3) * 

J 

+ \(° ijkl J o *^>/ ^'i)i 2 ("§;R) ' 
which with (|9.26|) precisely gives the strong convergence 

^ w in L 2 (0,T;H\n f ;R 3 )) (9.28a) 
iDi(t) -► t5(t) in L 2 (0;M 3 ) for any t G [0,T] , (9.28b) 



ft /■' 
/ wj_ 

Jo m ; 



I w in iJ^OgjR 3 ) for any i G [0,T] . (9.28c) 



Step 6. Uniqueness. 

Now, to prove uniqueness, let us assume that there exists another solution w' to (|9.1j) . 
such that w' G L 2 (0, T; Vo(-)), w£ G L 2 (0, T; L 2 (ft; M 3 )). By denoting foo = w - w' , we see 
that Sw G L 2 (0,T; V v ) is a solution of 

00 (5w t , <^)i,2(q. R 3) + K s fc<^V , a|<^>, s ) i 2 (n /. R3 ) + (c yW J Sw k , t , ^ )£2(ng ; R) = 0, 
(ii) «5w;(0) = in O . 

Since 5w(t-) G L 2 (0, T; V{;(-)), we can use <5w; as a test function in (i), which gives a.e. in 
(0,T) 

Id 
2~dt 



II^IIl2(Q;K3) + (c ijH / 5w k ,l , 5w\j ) L 2 {n s. t 
J 

+ u(a r k Sw, r , a s k 5w, s ) L2{n f. RS) = 



which, with the condition 8w(0) = 0, precisely proves that 5w = 0, establishing the unique- 
ness of such a solution. □ 

We will also need information on wtt- 

Theorem 9.2. Let w G L 2 (0, T; V«(-)) denote i/ie unique weak solution of \9.2T\) which 
is ensured to exist by Theorem 1.9. it TTien ^ G L 2 (0, T; iL^ (f2; R 3 )) and tD^ G W([0, T]). 
-2/n t. r?rnf-- 



Furthermore, q t G L 2 (0, T; L 2 (U J Q ; 
Proof. Step 1. Limit as e -> in (j9.17jh 
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In order to get an estimate independent of e from (|9.17|) . we integrate by parts in time 
to remove the second time derivative on q^' 

J 

+ q tt (t)[2(4) t w t l, j +(4) tt w e i , j ](t) 

- g x [2(5^(0)^ +(5?)tt(<W>i ], (9.29) 

1 f 

(g et (0) = q\ [ (a])j(0)uQ,j +a\ (0)w\,j ] = q\ in S7q by our compatibility conditions on 

the initial data) , from which we then infer by the regularity of a and of w e that 
1 /** 1 

^H^ttlli 2 ^;!^) + ^ y ifl r k w €tt , r , alw ttt , s ) L2 ^ n f. R ^ + -(c ljM w^a,w e l,j) L 2^s o . R) 



<cc 5 J o ku\\ 2 L2{U f. A) + s j Q llv^ll^ ( ^ ;R9) +5|| fe (t)||^ K;R) 

+ C C 5 ||V^(t)||2 2(n , ;R9) + (7 C 5 ||Vw £ (t)|£ 2(n /. R9) + C s CN{v ,f) 2 



where <5 > is arbitrary. 
Thus, for 5 small enough, 

\\ w ett(t)\\l2 {um + J HVWettll^a^/jjjB) + ll V ^ f (t)|li2 (ca . R 9) 



< C 



ll^lli 2 (^;R)+ 5 ll^(*)lli 2 (^;m) + (5 ^ll V ^(*)ll^(n^) 



+ CiY( Uo ,/)^ + C / |k e «||i 2(n . R3) . (9.30) 



By the Lagrange multiplier Lemma 18.61 we also have that 

|2 i ||V7„.. C+MI2 , /=»IIY7„.. 1 2 



+ d ll^)lli 2( n/ ;R ) + ll v ^(*)ll^(ng ;R9 ) + ^o,/) 2 ] , 
and thus with (|9.24|) . (|9.25|) and lj9.M()j) for a choice of 5 > small enough, 



+ WVMQWhwm + <5 ^(«o, f? + c J \\wett\\ 2 L 2 (n . m } 

which by Gronwall's inequality and (|9.24|) . gives 

^tf L2(Q f. R) < C[N(u J) 2 + [ \\w ett \\ 2 LHum } , 
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and thus, 

rt 

|2 /"* f A7" / „ , -f\2 i I|V7„,. ^+\l|2 i / ||„,, ||2 



MW L2(u f, R) < C [ N(u jy + \\Vw £t (t)\\ 2 LHQ f. A9) + J Q lhe«|li 2(n;R 3) ] • (9-31) 
We then infer from ()9.3(Jj) that 



K»(i)|| L 2 (n . R3) + J l|V^ £tt || L2(n / ;R9) + l|V«; £t (t)|| L 2 (n g. K 9 ) 

< C [ \\^ et (t)f L2{Q f ;M9) + N(u ,ff + f \W tt \\h {nv 



By the Gronwall inequality and (|9.24|) . we first get an estimate on j \\'Wett(t)\\'L2(Q-M, 3 ) w hich 
implies in turn that 



W ett(t)\\ L 2 (nm + / \\Vw ett \\ L Qf +\\Vw et (t)\\ L 2 {nSo 



< C \\V We M\l Hn>9) +CN(u jf. 

Step 2. e-independent estimate for w ett on [0, T]. 

In the same fashion as we derived (|9.24|) from (|9.22|) . we can deduce that for all t £ [0, T], 



IK tt (i)||i 2( fi ;R3) + J l|V^ e «||^ (n /. R9) + ||V™ et (i)||£ 2(ng . R9) < C N(u Q ,f) . (9.32) 

From (|9.31|) and (|9.32|) we then infer 

\\^(t)\\l 2(Q f ;M) < CN(u ,ff. (9.33) 

We thus deduce that for the choice e = — there is a subsequence, still denoted wj_, 
such that 



Wj_ - 

m l 


+ w in L 2 (0,r;F x (O;I 


t»)) 


(9.34a) 


W l - 

m l t 


±w t inZ/ 2 (0,T;ir 1 (n;] 


R 3 )) 


(9.34b) 


W 1 

m l tt 


*w u inL 2 (0,T;L 2 (O;: 


R 3 )) and in L 2 (0, T; JET 1 ^; M 3 )) 


(9.34c) 


qj- - 

m i 


+ q in L 2 (0,T;L 2 (^;I 




(9.34d) 


m i t 


inL 2 (0,T;L 2 (^;] 


R)) . 


(9.34e) 



Step 3. Initial condition for Wf 
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By the weak convergence in (|9.34|) . we infer from (|9.12|) that for each test function 
eL 2 (0,T;ftf(n;R 3 )), 

T r-T 

(W tt , <^)L2(n;K3) dt + V ({a r k a s k W, r )t, <t>, s ) L 2 (n /. M 3) dt 

f T f T 

+ ^ (cV kl w k , l , ( f>\ j ) L 2 (Q s m dt- J ({alq) t , <P\j) L 2 (n f m dt 

(^,0) L 2 (f7 / ;R 3) + £2 (n g. ffi 3) . (9.35) 



o 



Now for the initial condition, we notice that w t G C°([0, T]; L 2 (^; M 3 )). From the following 
identities which hold in L 2 (S7;IR 3 ), we find that 



w t (t) = w t (0) + / Wtt, w e (t) =Wl+ w ett . 

Jo Jo 

We deduce from the weak convergence of / w ett to / Wtt in L 2 (0, T; L 2 (f2; M. 3 )) that 

Jo Jo 

w t (0) = w x in L 2 (n-M. 3 ). 

Step 4. Strong convergence: the easy cases. 

We will also need the fact that the weak convergence in (|9.34|) is in fact strong. Notice 
that since wt S L 2 (0, T; Hq (S7; M 3 )), we can use wt in (|9.35f) to get for any t € [0, T], 



2ll^(*)lli 2 (n ; K3) + u J ((d r k a s k w, r ) t , w t , B ) L 2 [Q f. R3) - J ((a\q)t, ) L 2 (n /. a) 
+ ^(c ijkl w k :l (0)i 2 (ng;K) = jIHHi^nsR 8 ) + 2^ fe ^o>i ."oSj )i»(ng;R) 

+ / (^>^) L 2(q/ ;M 3) + (/ti*t)L 2 (ng;R3) • 







2 (^o ; R3 ) 

Since w(t) 6 Vg(i) in [0,T] implies (a^)^u; l ,j = —a\w\^, we then deduce 



1 ,2 



1 11 

+ ~{c l:>kl W k ,l (t),w\j (t)) L 2 (n s. R) = -|K||f2 (n . R 3) + -{ C l:>kl U k Qd ,U \j ) L 2 (n s. R) 

+ J (Ft, *i) L 2 (n /. R 3) + (ft, ^t)L2(ng ;M 3) ^. (9.36) 
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Similarly, 

^Kt(t)Hi 2( n;R3) + u J Q ((.°- r k d s k w € , r ) t , w € , s ) L 2 {Q f. R3) ~ {{a\<le)t, ™el,j) L 2 {Q f. R) 

+ -{(* ikl W e k ,i (t),W e \j (t)) L2( ng;M) = ^H^llll 2 ^ 3 ) + ^^^^l .Wy )L2(ng;R) 

+ / (^)^et) L 2( C /. M 3) + (/t)^et)L 2 (ng;R3) (it. 
J 

From the definition of q £ , (al)tw e l ,j = —a\w^j —e (q et — q\), and thus 
1 /"* 

+ "J ((a r k a s k ) t w e , r , w et , s ) L2{n f. m - J ((aj)tqe, w e \,j ) L 2 {n f. R) 
+ J Q ((al)tQeu <ij ) L 2 {n f, K) + (Qet, Qet ~ Ql) L 2 {n f. R) 

+ -(<? jkl W e k ,l (t), V)j,j (t)) L 2 (n s. m = - HwiHia^RS) + -{c ijM 4, l ,U i , j )L2 (ng;ffi) 

+ Sq ^^LH&om + f Q (ft> W et)L2(n°;R3) ■ (9-37) 

By integration by parts, since q e (0) = qo, 

J ({ajUe, W e \,j ) i 2 (Q / ;M) dt = ~ j ^Me)t, ^e\j) L2{n f. R) dt 

+ ((al)tqe{t), We\j{t)) LHn f m 

- ((aj)t(0)q O , «0 < ,i) L2(n /. R) . (9-38) 

Now, since g e (t) = qo + qet dt we deduce that for any i G [0, T], q e (t) ^ qo + Qt dt 

Jo Jo 

rt 

in L 2 (Q£;IR), which proves that q(t) = q + / <? t , and thus that q £ C°([0,T]; L 2 (%;R)), 

Jo 

with q(0) = qo- Furthermore, for any t G [0, T], 

q e (t) ->> q{t) in L 2 (f^; M) as e -»• . (9.39) 
Similarly, since w e (0) = w(0) = uq we have that for any t 6 [0, T], 

w e (t) -» w(t) in i^^M 3 ) as e -► . 

Moreover from 
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we infer from the strong convergence in ()9.28[) and the weak convergence in (J9.34|) that for 
any te [0,T], 

/"* 1 

from which we obtain the strong convergence 

wJt) -> u)(t) in iJ^O^M 3 ) as e = — -> . (9.40) 

mi 

Thus, from Q9.38|) . the strong convergence in (|9.28|) and (|9.40j) together with the weak 
convergence in (|9.34|) and (|9.39j) shows that 



rt rt 

J Q •• -■ ~ -"-\"am j Q " -' -■--■^v.«ow mi 



((2?) t g 6 , Wel,j) L 2 [n f m I ((3?) t g, ^ ) L2(n /. R) d* as e = -> . (9.41) 



From (|9.41[) . the weak convergence in H9.34|) and the strong convergence in (|9.28|) . we 
then deduce from (|9,36|) and ()9.37j) . that as e = jjjj- — ► 0, for any t G [0, T], 

1 /"* 

2lKt(i)IU 2 (n;R3) + 17 y K a fe w et,r, ^et,s ) Z 2 (n /. R 3 ) ^ 

+ I( c ^^,,(t), «,* -(t)W ; 



l /"* 

1 

2 l 



+ -(c«*'ui*, I (i), 



which implies the strong convergences 

to i (t) -> u>(t) in ^(tyM 3 ) for any t G [0,T], (9.42a) 

in L 2 (0,T; H 1 (Q f -M. 3 )) (9.42b) 

to i (t) -> w t (i) in L 2 (ft;M 3 ) for any f € [0,T] . (9.42c) 

From the strong convergence in (|9.42j) and Lemma l8.6[ we also deduce that 

II* " 9llza(0,r;La(^;R)) 0, as e ^ . (9.43) 

Step 5. Strong convergence: the more delicate case of w t t- 

Our main difficulty results from the fact that we cannot directly obtain an energy inequal- 
ity for wu (from the limiting weak form of the twice time-differentiated problem). Rather, 
our starting point will be (|9.17|) . from which we will get by weak lower semi-continuity the 
desired inequality, provided that we can prove that w ett — > Wu in L 2 (0, T; L 2 (Qq ; M 3 )). To 
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prove this result, let us first remind the reader that w e n satisfies for all <j) G L 2 (0, T; Hi (fi 



+ 



T 



l Wet,i , ( l ) \j)L 2 (nf ) ;R) dt - / ((a k q e )u, (f> jl) L 2^ n f. u ^dt 

J 

T 

(F it ,0) L2(n /. R3) + ^)i^(ng;R3) . (9.44) 
From the bounds associated to the weak convergence in (|9.34|) . we then see that 

Kttt(*)llv 6 (t)' dt <C, (9.45) 



where C denotes a constant which depends on the data, the smoothing parameter implicit 
in a, but not on the penalization parameter e. In the following, this letter will denote a 
generic constant depending on these variables. Let us fix 5 > and let 

Of = {x G dist(x,r ) > 5, dist(x,dO) > J} . 

Let us then denote on [0, T], 

n,(t) = • 

For each i G (0, T), we have the existence of 5i > such that 

W e(t-st,t + st), n 2 s(t) c n s (t'). 

By a simple change of variables and (|9.45|) . we then get 

||deto«> ett t o r^\ L 2 (t _ Stit+St . H i diJ nMt)m') ~ G ■ (9 - 46) 

We set 



u e = det a w ett ojj 1 



From (|9.46|) and (|9.34|> . we then get 



l M ^llL2(t-5t,t+<5t;^ div (n 2i (t);R3)') < C 
ll w e|ll,2^- ( Jt,t+ ( K;^(n2«(t);7 n ' 3 ^ - ^ 



Thus, from the classical compactness results (since the domain Cl2s(t) is fixed on (i — 6t, t + 
fit)), 

u e -> u = detaw i4 o fj- 1 in L 2 (t - 6t,t + 5t; L 2 (n 2 s(t);R 3 )) , 
which obviously gives (since on (t - St,t + <5i), ^(i', n 2( j(i)) C 

w ttt -► «fc in L 2 (t - <5i,i + ft; L 2 (o£; M 3 )) . 

By a finite covering argument, and the L°° bound 1)9.32(1 of w ett in L 2 (Q; M 3 ), we then infer 
that 



limsup / KttlljW.RSj d *= / * • ( 9 - 47 ) 
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Successively from the Cauchy-Schwarz and Sobolev inequalities, 

< c \n f nnf \ N(u ,f) 2 . (9.48) 

Prom (|9~T7j) and (|9.48|) . we then infer 
f T f T 

hmsup ^ Kttlli a(n /. R a ) *<y o ll%& (n /. R3) di + c|n£nnf| iVN,/) 2 , 

which immediately shows that 

fT 
'0 



hmsup / IK«ll^ (n /. R3) ^< / IKII^n^s) dt 



thus establishing the strong convergence as e = — ► 0, 

-> % in L 2 (0, T; L 2 (^; M 3 )) . (9.49) 

m ! tt 

Next, we restrict our test function (p to be in the space {(ft E Vy(t) | = on f^}. For all 
such test functions and for a.e t £ (0,T), G Vs(f), 

(%(t) - W e «(i),0)L 2 (n ; R3) + K(Ofcaf(* - W e , r )) t (t), JS ) X 2 (n /. R 3) 

- ((Oi)t(g - g 6 )(t), ) L 2 (n /;M) = ; 

thus, the second Lagrange multiplier Lemma 18.81 ensures us, from the strong convergence 
in (I9~i31) and (l9~i2cl .that 

1 f 

where qt = qt y / qt det V77, and a similar definition for q e . In the following, we will 

If If 

denote c = — — / qt det Vt? and c e = — j- / q et det Vfj. 

Step 6. An inequality for wu with a constant independent of the mollification 
parameter. 

Now, from the weak convergence (|9.34|) and the compactness of the trace operator, we 
then infer that as e = — > 0, 

IK* - w e)tt\\ 2 L 2 (0tT . L 2 (To . m) -» . (9.51) 
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We now note that from (j9~TT|) and for any < t < T, and < St < Mm(t, T-t), 



t+st 



2 Jt- 
1 



+ 



-St 



\ W ett\\L 2 (fl;) 



+ V 



t+St nt' 
t-St Jo 



t+st 



2 J t-St 

t+st 



[c ijkl w et h \i ,w et < i ,j) L 2 (ng . I 



{0.k w ettir > &k w £ttis )i2(n/.R3) 
i+<5t /-t' 



t— 5* ./o Jn£ 



Qet [2 (a^t^et^+Ca^ttw/^ 



g et [2(5^)^4, j +(aj)ttw € \ j ] t 
t+st r t' 

(al)tqetWelt,j 



t-St Jo Jn? 



+ 2 v 



t-st Jo Ju f 

t+St rt' 

t-St Jo 



'M)ttqeWel t ,j +V 



t+St ft' 



t-St JO 



((a r k a s k ) tt w e:r , Wett,s) L 2 (n /.j 



)t w etir j w tttis )x,2(q/ ; k3) 
«t+<5t /•*' 



nt+Ot rt rt+St rt' 

<C5tN(u ,f) 2 + / (Ftt,w ett ) 2( f . u3) + / {ftt,Wett)iJ>Ki°-m- (9-52) 

Jt-<5t JO JO 

The first three terms of the left-hand side of this inequality will be dealt with by weak 
lower semi-continuity. By the weak convergence in (|9.34|) and the strong convergence in 
Q9.28JI . ()9.42cj) . q9.4U|) . and ()9.43|) . we infer that all of the remaining terms, other than the 

term Jq J n f (ai) t q et w t l tt ,j , converge as e — > to the same expressions with the limits 

w and q replacing w e and q t . From the definitions of c and c e , we have that 



t+St ft 



t-&t Jo Jo, f 



(al)tqe t WeU,j 



t+St rt' 
t-St JO 



Sal)t{qe)tw e l t ,j +c € I f (al)tWett>j 



From the strong convergence (|9.5Uf) and the weak convergence Q9.34J1 . we then deduce that 
the first term of the right-hand side of this inequality converges as e = — — > to the 
corresponding term where qt replaces q €t and wtt replaces w ett . 

For the second term of this right-hand side, we notice from a spatial integration by parts 
(since c e depends only on the time variable) that 



t+St ft' 



t-St Jo 



,a J ihu>ett,j 



rt+St rt' 




/ / Ce 


In 


It-st Jo 





and thus from the weak convergence in (|9.34j) and the strong convergence in ()9.49[) and 
(J9.51|) we then get the convergence as e = ^- — > to the corresponding term where c 



replaces c e and ibu replaces w ett . This implies that as e = ^- — > 0, 



t+st ft' 



t-st Jo Jn{ 



{al)tqe t Welt,j 



t+St ft' 



t-St Jo Jn f 



(dl)tqtw l tt ,j . 



Consequently, all the terms, except the three first ones, appearing in the inequality (|9,52|) are 
convergent < 
by q and w. 



convergent as e = ^- —>■ to the same expression, where q e and w e are replaced respectively 
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By weak lower semi-continuity for the first three integrals, we then infer that as e = 
^- — > the same inequality as the previous one holds with w and q replacing respectively 
w e and q e . By dividing those integrals by 25t and passing to the limit as 5t — > (which is 
possible a.e. in (0, T)), we then get that a.e. in (0, T), 

t 



T, \\Wtt 



2 n--\-/nL2( n . R 3) + y (a r k Wtt,r, ^k^tt,s) L 2^ n f. 



+ \{c i3kl w k t ,l (t),w l t ,j (t)) L 2 {nsm - [ f q t ^(aDtwlj +{aP i ) tt w\ j ] t 
z Jo Jn' 



&(t) [2(a?) t T5j )i +(3?) tt tS t , i ](<)-2 / / fo?)*^,,-- / / ral)ttqwlt,j 



u JO Jf 1 ./n ./<)' 



" y ((5fc5fc)tt"V, *tt 5 s) L 2 (n /;R3) + 2 " y ((a r k a s k )tw t , r , w a ,s) L 2 {n f 

rt r t 

2 



<CN(u ,f) + / (F tt ,%) i2( j J / ;R3) + / (/«,t3tt)z^(ngiR3)i ( 9 - 53 ) 

>/ " 

where (we recall) C does not depend on the smoothing parameter of a. □ 

9.5. Regularity for w and its first and second time derivatives, dependent on the 
regularization parameter of a. As discussed in the introduction, we shall focus on the 
regularity near the interface, which will provide us with the trace estimates on the interface. 
Elliptic regularity for the Dirichlet problems will then yield the full regularity result in each 
interior component. In this subsection, C continues to denote a generic constant which 
depends on the same variables as C and C(M), and additionally on the regularization 
parameter. In Section I1U1 we obtain estimates independent of n, by interpolation mainly, 
which requires us to know a priori that the solution is smooth (without using the estimates 
that we get in this subsection, since they blow up with the regularization parameter). 

Recall that we have already shown that w £ L 2 (0, T; V«(-)), w t € L 2 (0, T; H^(Q; R 3 )) and 
w t t 6 W([0,T]), and that both q and q t are in L 2 (0,T; L 2 (n f Q ;R)). 

The missing regularity results will be recovered using difference quotients. Recall that if 
we consider the partition of the space M 3 formed by the two half-spaces := {(x 1 , x 2 , x 3 ) E 
M. 3 | x 3 > 0} and M 3 := {(x 1 , x 2 , x 3 ) £l 3 | x 3 < 0} and the horizontal plane {x 3 = 0} with 
the usual orthonormal basis (ei, e2, 63), then we have 

Definition 9.2. The first-order difference quotient of a function u of size h at x is given 
by 

u(x + h)- u(x) 
D h u{x) = , 

where h is any vector orthogonal to e%. The second-order difference quotient of u of size h 
is defined as D_hDh,u{x), given explicitly by 

_ _ , N u(x + h) + u(x - h) - 2u(x) 

D_ h D h u(x) = —2 . 

\h\ 2 

We will denote 

V Q u = (n,i , u, 2 )• 
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Letting A l i = a\a d k , we write the weak form as 

(W t , ^)L2 (n . M 3) + u(A v W,i , <f>,j ) L2(n /. M 3) + (C 1 ^ 1 [ W k ,i , 4>\j )L2 (Q g. K) 

J 

- (9.4^ fc »i)L2 (n /;R) = ( F ^)l2(^;K3) + (/,<A)L2(ng;R3) 

for all G ^(fijM 3 ), for a.e. < t < T. 

Next, assume that ft = 5(0,1), the unit ball centered at 0, and that = {x G 
5(0,1) | x 3 > 0} and fig = {x G 5(0,1) | x 3 < 0}. Select a smooth cutoff function ( 
satisfying 

1 q 

C=l on 5(0,-), ( = on M 3 -5(0,1), and < C < 1. 

Let <p = D_ h (( 2 D h w); then clearly <j) G 5,j(fi;IR 3 ) for a.e. t G [0,T]. We may thus 
substitute <\> into the above weak form to obtain 

A 1 + A 2 + A 3 -A 4 = 5, 

where 

A 1 = (D h w t ,C 2 D h w) L 2 {n . R3) , 

A 2 = u{D h (A^W,i ), (C 2 A^),j ) L 2 {Q f. R3 y 

A 3 = (D h (c^ hl w^O^^DhW^j)^.^, 

A A = (q,a l k (D^ h [C 2 D h w k ]), l ) L2{n f m , 
B = (F,D- h (( 2 D h w)) L2{Q f. R3) + (D h f,( 2 D h w) L 2 (n s. R3) . 
For the first two terms, we easily find that 



A 2 > C\\CD h Vw\\ 2 ^ - C\\Vwf 



For the remaining terms, we shall use the notation C h (x) to denote C(x + h). Whereas 
the coefficients of the elasticity tensor are constant, one should keep in mind that our 
assumption on the domain comes in fact from a change of variables which produces a non 
constant elasticity tensor. It is the integral below with the D^C term which necessitates 
the hyperbolic scaling of our functional framework. 

Expanding A3, we have that 

A3 = (C 2 c h : D h J Vw,-D/ l Vw) £ 2( n g ;R 9) + (C 2 D h c : Vw, D h Vw) L 2 (n s. K9) 
+ (2(V(®D h w,c h :D h ! Vw) L2(n g. M9) + (2CC,j D h c ijkl f w k , t , iW^ngjR)- 

J J 
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The second term on the right-hand-side is 

{D_ h (( 2 D h C: I Vw),Vw) L 2 {n s. R9) 

J 

= (l( 2 D h C]- h : D-h f Vw + D_ h {( 2 D h C) : [\w,Vw) L 2 {n s R9) . 
Jo Jo 

Hence for > 0, we see that the A3 term yields the following inequality: 
\A 3 - \j t (C 2 C h : jf D h Vw,J* D h Vw) L 2 {n s m \ 

< e ^ D - h J Q v *lli 2 (n§;R 9 ) + c o\\Vw\\ 2 L 2( n ..j i ? ) + C\\ J Vt2)|||2 ( n Si R9). 
For the A4 term, we have that 

M = (q , a{[( 2 ]- h D_ h D h w k , l +a{(D_ h ( 2 )D h w k +[2CC,j }- h D_ h D h w k a[ 
+ 2a[D- h (tt, l )D h w k ) L2{n f ;R) . 

By the divergence-free condition, w G V{i([0,T]), we get in S7q 

= D_ h {[a{} h D h w k , l +D h a l k w k , l } 

= ~a l k D„ h D h w k :l +D_ h [a{] h D h w k , l +[D h ~a{]- h D_ h w k , l +D_ h D h a l k w k ,t , 

allowing us to eliminate the first term appearing in the expression of A4, which gives for 
6 > 0, 

\M < e\\CD h u, k , \\ 2 L2(nlm + c e \\q\\ 2 L2{n! . R) + c\\Vw\\ 2 L2 ^ m . 

Finally, 



\B\ < 0[ ||V^||i 2(ng;ffi9) + IIC^V*||^ K;R9) ] 

+ C \™% {<M9) + C 6 [ Wffn^m + 11/11^3) ]■ 
Choosing > sufficiently small, we have the inequality 

j t (\\(D h w\\ 2 L2{n . R3) + (C 2 C h : J* D h VwJ* D h Vw)„ in . m ^ + IIC^V^^,.^ 

< C fa f ^ V *IIl 2 (^;R9) + II J V *llL2(ng;R9) + 1 1 V ^ 1 1 ^ (f^g ; R9) 

+II v *IIl 2 (^;r 9) + W^Unf-R) + WK^nf-RS) + WfWnH^m) ■ 

From Gronwall's inequality, it follows that d a djW G L 2 (0, T; L 2 (W ; R 3 )) where V-'' = {2: G 
B(0,i) I x 3 > 0}, and where a = 1,2 and j = 1,2,3. Hence, € L 2 (0, T; ^(W; M 3 )), 
so that by the trace theorem we obtain d a w G L 2 (0, T; i?°- 5 (V r/ n {x 3 = 0}; R 3 )). Thus, 

u) G L 2 (0,T;H 15 (V f n {x 3 = 0};R 3 )). (9.54) 
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/' 

J 



(with an estimate which blows-up as the mollification parameter n — ► oo). Similarly, 

ft 

w £ L°°(0, T; 5 L5 (V S n {x 3 = 0}; R 3 )) , (9.55) 

/o 

where F s = {x G 5(0, ±) | x 3 < 0}. 

We now drop the assumption that f2 is the unit ball, and once again assume it is an open 
bounded subset of R 3 with all of the smoothness assumption stated previously. We choose 
any point xq £ Tr, and assume that 

n B(x , r) = {x G B(x Q ,r) \ x 3 > ^(x 1 , x 2 )} 
Q s n B(x , r) = {x £ B(x , r) \ x 3 < ^(x 1 , x 2 )} 



for some r > and some smooth function 7 : R 2 — ► R. We define the following change of 
variables: 



y * = X i =: &(x), i = l,2 

y 3 = x 3 - 7(x\x 2 ) = $ 3 (x), 

y = ®(x). 

x* = y* =: i = l,2 

x 3 = y 3 + 7 (y 1 ,y 2 ) = * 3 (y), 



and write 
Similarly, we set 

and write 

x = #(y). 

Then $ = and the mapping x 1— > <fr(x) = y straightens out Tq near xo, and det $ = 

det^ = 1. 

We can assume = $(xo). Choose s > so small that 5(0, s) C $(B(xo,r)). Let 
= «>(*, g / (t,y) = ^,$(y)), f'(t, y) = f(t, d» (y)), 

Then it/ and io{ are in L 2 (0, T; 5 1 (r2 / ; R 3 )), where Q' = 5(0, s). We also set 

n{' = 5(o, s) n {y 3 > 0}, n s ' = 5(0, a ) n {y 3 < 0}, 

Then, since (w, q) satisfy the weak formulation, applying the change of variables, we see 
that (w',q') satisfy 

r-t 



H, cf)') L 2^, R3} + v{a' ij w',i , <j>',j ) L2(n f' m + (c Hjkl w' k ,i , 0*j ) L 2 {Q g/ ;R 3 ) 

- W, [a[ o W*, r ) L 2 {n f> m = (/', ^)L»(n';R3) (9-56) 

for all <// G 5 1 (O';R 3 ), for a.e. < t < T, where 

a' kl =A^o^ g k g), ^rks = c ijkl g s g r ^ g^y) = ^y)]~\ 
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It is easy to verify that both a' and c' retain the uniform ellipticity conditions of the 
original operators A and C; Moreover, w' satisfies the divergence condition a\ o $>w h ,j = 
in [0, T] x f2'. Thus we may apply the results obtained above for the case that the domain 
is the unit ball to find that 

w' G L 2 (0, T; H x *(y f ' n {x 3 = 0}; M 3 )), 

/ w' G L°°(0, T; H i r °(V s ' n {x 3 = 0}; M 3 )), 
J o 

where V f> = {x £ B(0, f ) | x 3 > 0} and V s ' = {x G B(0, § ) | x 3 < 0}. Consequently, 

w G L 2 (0, T; H 1,5 (dVf n To;lR 3 )), f w G L°°(0, T; F L5 (dy s n T ; K 3 )), (9.57) 

jo 

where V* = ^(V f> ) and V s = ^(y s/ ). 

Since To is compact, we can as usual cover To with finitely many sets of the type used 
above. Summing the resulting estimates, we find that we have for the trace on To 

w G L 2 (0, T; H 15 (r ; K 3 )), (9.58a) 

f w G L°° (0,T; H 1,5 (Tq;M. 3 )). (9.58b) 
Jo 

Converting the fluid equations into Eulerian variables by composing with 77 , we obtain 
a Stokes problem in the domain ^(Oq): 

— uAu + Vp = f — wt o f]~ x + vaj,j or]^ 1 u,i ~P (5f ),j °V~ : (9.59a) 
divu = 0, (9.59b) 

with the boundary conditions that u = 0on fj(dQ) and that u G L 2 (0, T; H 1 - 5 (fj(T ); M 3 )), 
where u = wofj -1 andp = qofj" 1 . Since the domain 7?(Qq) is of class H 3 , by the elliptic regu- 
larity of Q2], (l93HaT) implies that u G L 2 (0, T; H 2 (fj(n f ); R 3 )) andp G L 2 (0, T; fl^fij);!*)). 
It follows that 

«ieI 2 (0,T;5 2 (^;R 3 )), g G £ 2 (0, T; H 1 ^^; R)). (9.60) 
Similarly, elliptic regularity of the elasticity problem shows that 

/ wa°°(0,T;ff 2 (liJ;R 3 )). (9.61) 

Next, we consider the weak form for the time derivate wt for all G Hq (O; R 3 ): 

(*«,<A)i2(C;R3) +^([^ rS *,r]t,0, S ) L 2 (Q / ;R 3) + (c Vkl W k ,l , )i2(ng;R) 
- ([^]t'^j)L2 (Q / ;R) = ("^'0)^2(0/^3) + (/t, 0)^3(05^8) ffl - e - ^ G ( ' T )- 

Expanding the time derivative, we see that there are two additional terms in the weak form 
given by (A r t s w, r , cp, s ) i2(n /. R 3) and {{at)\q, ) L 2^ n f. R3 y These additional terms are easy 
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to handle, and by letting (p = D-h{( 2 DhWt), and following the identical procedure as above, 
since we also already know that wtt G L°°(0, T; L 2 (S7; M 3 )), we find that 

w t G L 2 (0,T;H 2 (n f ;R 3 )), q t G L 2 (0, T; H 1 ^^; R)), * G L°°(0, T; F 2 (^; M 3 )). (9.62) 

Because of the assumptions on the forcing and these estimates for -tit, we may improve 
the regularity results H9.6UJI and (|9,61|) , We apply the identical procedure, but this time we 
use 4> = D_i l D) l (( 2 D_i l Di l w) as the test function. We find that 

w G L 2 {0,T;H 3 (n f ;R 3 )), q G L 2 (0,T; H 2 (Q f Q ;R)), [ w G L°° (0,T; H 3 (Qq;M. 3 )). (9.63) 

J o 

Moreover, [|(i2j,g)[|^ r < CN(uo, f), where the constant C — > oo as the mollification param- 
eter n — ► oo. 

In the following section, we will use a different form of (|9.56j) . If we denote by £ a 
smooth cut-off function, equal to 1 in a neighborhood of contained in $7' and outside 
0', and denote W = ( 2 w', Q = ( 2 q', bj = af o $ ^, C 4 -^ = c /yfc ', we then obtain for any 

G i7 1 (M 3 ; M 3 ), 

(W t , ^)i2(M3;K3) + ^ (trffiWik , <P,j )l 2 (R1;M3) + / , y )l 2 ( R 3_ ;R 3) 

•/ 

(Q> &fc^ >r )l 2 (R3 ;R) 

= (-^1) ^)l 2 (R3 ;R3) + (Hj, <f,j ) L 2( R 3 R 3) + (-^2, V)l 2 (R 3 ;R3) + (-^j ! )l 2 ( R 3 R 3) , (9.64) 



where 



= C 2 ^ - * C 2 ,i w H , k +q'b%C 2 ,r , (9.65a) 



'J "i 



Hj = v bjb k ( 2 , k w', (9.65b) 



c 2 f ,i_ c mi C 2 /'„/*,, (9.65c) 



K) = C ljkl C 2 d j l w' k . (9.65d) 

Moreover, satisfies the divergence condition 

j- = a = ^C 2 ,j M 7 * m [0, T]xi 3 . (9.66) 

Note that we consider the above inner-products over all of M 3 since W and its derivatives 
are compactly supported in 0'; the contribution outside of Q! is zero regardless of the way 
in which we extend b and g to [£l'] c - This same remark also applies to (j9.66|) . 

10. Estimate for (J7TTJ): the case of the actual coefficients 

10.1. Energy estimate for wu independent of the regularization parameter for a. 

We are now going to use the regularity results (|9.62|) and (|9.63p in the energy inequality 
(|9.53|) (which was bounded by a constant that does not dependent on the mollification 
parameter). Our approach now will be to use interpolation inequalities to obtain an estimate 
which is independent of the regularization parameter. 
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This section will be divided into eight steps, each of which is devoted to the estimation 
of the various integral terms in (|9.53j) . 

In what follows, 5 > is a given positive number; the choice of S will be made precise 
later, as it will have to be chosen sufficiently small. 

Step 1. Let I\ = I I qt(al)tw\ t ,j . Then, by the Cauchy-Schwarz inequality and by 
Jo Jtt{, 

interpolation, 

h - S J \\ V ^Unim +C5 l ll 5t H^(^ ;R9 )ll^ll^ (n / ;]R) 

where we have used (|6.3j) for the L°° control of at in H . Thus, 

ft 

id./) " Jo 

By Lemma 18.61 applied to (|9.35l) , and 1)6.3(1 , 

l^&(^;R) ^ C [ H^Hi 2 (n;R3) + II? «tlli 2(n /. R9) + ll^ll^ (n / 

+ ll/tllL»(ng^3) + l|wt&i (n / ;m 3) + H*II^K;iR3) + IMI Wg;^) ] • (10.1) 

Thus, with wt(t) = w\ + wu, w(t) = uq + / wt and q = qo + q~t respectively 

Jo Jo Jo 

^(nfjR 3 ), # 2 (f^;IR 3 ), and fT^n^R), 



h < 8^ ||V^|| 2 L2(n / ;R9) + CsC(M)*V Utility Ijf H&lSn^)] 5 Ti • 
mma 18.61 applied to (|9.35l) . and (|6.3j) . 

II 2 s n r ILr, l|2 i ||^ ~ 1 1 2 , II p 1 1 2 

|2 , II ~ 1 1 2 i || ~.||2 



111 







+ CsC(M) T\ [ N(u , ff + T £ Utf HHn f m + [ »* 112 



II 2 



H 2 (Q^;R 3 ) J 



+ SUP 11*11^1(^3) + SUp ||%|| X 2 (nR 3) ]. (10.2) 
[0,T] V ' [0,T] 



Step 2. Let J 2 = / / ft(5?) tt u;J,j. Then, 



I 2 < / IK^)tt^.ill i f ( „/. R) llftllL«(n^^ 



o 
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Thus, using (|6.3|) for the L°° control of au in L 2 , 

h < 



which with wt(t) = w\ + I wu gives 

J o 



t 

2 



12 * Vo ^HHntm 

+ C 5 T\{ N(u , ff + f l|V^||^ (n , ;R9) +T J* ||V^f i2{n / ;R9) ] • (10-3) 




Step 3. Let h= / qtifll)tttw\j ■ Then, 



'o ./n£ 



2 IIV7„r.l|2 /o / ||= ||2 



/ 3 <*/ 11^11^(^^)11^11^;^) + ^/ II^H^K;! 







2 



ll ^ ll L 2( n/; R9) [ Ar ^'/) 2+T / II^H^K^)] 
+ C S J q \\qt\\ I | (n / ;R) 1 1 © 1 1 ^ (fl / ;R) , 

where we have used w = uq + / vbf Thus, by (|6.3[) . and Q1U.1JI . 

J 3 < «7(Af) [ N(u jf + T j* \\w t \\ 2 HHn f m ]+C s T\ ^ ||^|| 2 Hl(n / ;M3) 

+ ri [ n( U0 , f f +T f Utf m{nlM) + tJ* ||%||| 1(n / ;R3) 



+ T l H^ll^^^) + sup H^H^i^g.KS) + sup H^lli^^RS) ]. (10.4) 



t 

3" 



Step 4. Let h = I I q{al)uw\t,j ■ Then, 
Jo Jnl 



t rt 

2 i r< / lis l|2 ||~||2 



2 , /o -o/7i^\ / II ~i 1 0.5 n ~||1. 5 



- 6 1 llv^lli^ + ^Wjf ll<^ ;M) lk1^ 1(n / ; 



2 



- 5 / ll^lll,^. 



W (M) Ti [ JV(«o,/) a +T / lift 11^^, + f M\\ 2 HHn , R) ] • (10.5) 
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The next two steps will require the introduction of 8\ > which is different than 5 and 
will also be made precise later. 

Step 5. Let I5 = — qt{t){al) tt w l ,j (t). We first notice that 

h = ~ / f qt^ihiw^j (t) - ulj ) - I qtWiaDttulj . 

For the second term of the right-hand side of this equality, 

|2 , ri \\7. H„, I|2 





and thus by (|6.8|) . since T < Tm, 

- f^ f m$i)tt<i < SiUt(t)\\ 2 L 2 ia f m + C 5l C(M) N(uo, f) 2 - (10.6) 
For the other term, 



q t (t)(ai)tt(w\j (t) - ulj ) < 5 llftC*)!!^^ (W.7) 



+ ^l|5«(*)lli3 (n / ;R9) HV*- Vu \\l 6{nfo;R9) , (10.8) 
and thus, by the L°° control in L 3 provided by (|6.10|) . 

QtmaiMw^j (t) - uij ) < 5 \\qt(t)f L 2 {n f. R) 

+ C 5 C(M)T j\\Vw t \\l 1{nlm . (10.9) 
By (UnH, (tiTTKl) and (fTu~31) . we finally have 

h < (6 + [ N(u , ff + C(M)T [jf \\q t f Hl{Qim + jf ll^lli 1(n / ;R3) ] 

+ C 5 C(M)T l|V^||^ 1(o/ . R9) + C^C AT(^ , /) 2 ■ (10.10) 

Remark 20. -/Vote i/iai L 3 and L e in \1U. 8\) are limit cases for both \6.1U\) and the Sobolev 
embeddings in dimension three. In dimension > A, this would no longer be possible and we 
would be required to introduce a smoother functional framework. 

Step 6. Let Iq = — qt{t){a{)tw\-,j (t). Similarly to our previous step, we first notice 
that 



n-^ J of 

1 'o 
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For the second term of the right-hand side of this equality, 

f zi\.fri\„7,i . (+\ S A. Np;.cnM|2 iu.r+M|2 , n. II V7»r,. m l|2 



o 

and thus, 



^tW^'MOH.W < C*i iV(n ,/) 2 IIMII 2 ^,. 



+ ^ 1 [||v* 1 || 2 2K;M9) + r /Wll^/*., ] ■ (io.il) 



For the other term, 

'7,i\.i+\ _ fsiwrm-r,* . m * h^v-aii 2 



"0 



ft(f)((5|) t (t) - (^)t(O))^,, (t) < 5 ||ft(t)[li, (n / s 
and by ([577]) . 



+ Qll«t(i)-^(0)|li 8(n / ;R9) ||V^|| 2 i3(n , ;R9) , 



, ft(*)((Sj)t(*) " (^(0)X- (*) < * llft(t)lli 2(n /. R) 

+ C 5 C(M)r||V^(t)|| 2 3(no/;M9) 
In the same fashion as we proved ()6.1U|) . we use the L°° control in L 3 : 

n v ^(*)iii3 (n / ;K9) < iiv^iii3 (n / ;R9) + c [f Q iiv^ii^ 1(Co/;R9) + iiv^n 2 i2(Q , ;a9)J) 

which combined with the previous inequality provides us with 

ft(*)((2ft*(t) - (4)tm4,j (t) 



\l 



z 'IIMII^ + ^WrilVtSxii^ 



+ CsC(M)T[J q \\Vw t f Hl(n f m +l HVt5«|£ a(n , ;R9) ] • (10.12) 
By (|10.11j) and pU.12j> . we finally have that 



h< + [CN(u J) 2 + C(M)T / \\wtt\\ mm , 



Hlffl? ;R 3 ) 

+ c(m)t[ y o ii&ii^^ + ^ IKII^ (n , ;R3) ] + sup IKIIW) 



+ SUp ||^||jfl(f2g;K3) 



[0,t] 

t 



+ C 5 C(M) T [ N(u Q , f Y + / llVwttHj,^ + / \\^t\\ l m{ai . m \ • (10-13) 



t 

2 
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Step 7. Let I 7 



J ((a r k a s k )ttw, r , Wtt,s) L 2 in f. M3 y Then, 



h < sJ q Iiv^||^ (n /. M9) + c s J q ll(«fc«Dttlli 2(n /. R) llv*||^ li4(n /. R9) 



where we have used (|6.3|) for the L°° control of Ou, a>t and 5 respectively in L 2 , H 1 and H 2 . 
Thus, 



i 

2 



l7<$ I ||Vw tt || rW 



L 2 (n^;M9) 

t 



+ W(M)Tl [iV(no,/) 2 + T jf llVtOtll^^+jT IIV^^ ] . (10.14) 
Step 8. Let J 8 = - J ((a r k a s k ) t Wt,r , *tt, s ) L2(n /. R3) - Then, 

Is < 5^ \\Vwtt\\l 2(Q f m + C s ^ ll(^5|) t ||2 4(no/;R) ||V^||^ 4(n ,. K9) 



< tf* \\Vw tt \\ 2 L2(n>9) + CsC(M) £ l|V^|l° 2 5 (c / ;K9) l|V^|l^ K; 



Consequently, 



h<5 I \\Vw tt \ 12 



o 



+ C 5 C(M) T\ [NM+tJ q l|V^ t f L2(n / ;R9) +^ l|V^||^ 1(n / iR9) ] • (10.15) 
Step 9. Thus, from (|9.5^|) . and estimates (|1().2|) - (|1().15|) . we then obtain the inequality 

sup|K(i)|| 2 2(n . M3) + ^ 11^61^3) + sup ||t5t|||ri(n 5; R3) 



[0,T] JO ^ °' ; [0,T] 

< [C<5(l + C(M))iV(no,/) 2 + ^i(l+TC7(M) + iV( U o,/) 2 ) ] ||(«>,g)||| T 

+ Ci 1 JV(« , /) 2 + C 5l T3 ||(^,g)||| T + C 5 C(M)T3 (||(^,g)||| T + iV (u , /) 2 ). (10.16) 



10.2. Estimate of w t independent of the regularization of a. 

In this subsection, let Vl/j be one of the H s charts defining a neighborhood of Qq and let 
Wi = Ci w o Since the estimates that follow do not depend on the choice of we will 
denote Wi simply by W. 
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Recall that for all <p £ L 2 (0, T; 

J (W u ,<f>)L2(RZ;R 3 ) dt + V j Q (& r k b' k W, r )t, <^,5 )l2 (r 3_ ;R 3) dt 

+ j\c irks W k , r , <j>\ a ) L2(ffi3 _ m dt - <f>\j )LHRlm dt 



/ ( F U, 0)l2( K 3_. M 3) + (H it , 4>,i )l 2 (R^ ;R 3) & 
J 

+ / (-^2t, <A)l2(r3 ;K 3) + (-^it,0,i)L 2 (R 3 _;K a ) ^ ■ 
■/ 



With the choice of = D-hDhWt in this variational formulation, which is possible since 
w t € L 2 (0,T;H^(n-,R 3 )), we then get 



i||^m(T)||| 2(R 3. R3) + ^^ (6^D ft W t , r , ^^ t , s ) L2{R 3. M3) 
+ ±(C irks D h W k , r (T), £> h W^ (T)) L2(K 3 ;R) - £(D h (b(Q) t , D h Wt,j ) l2(r3+;R) 
+ uj\D h (b r k b s k )W t , h r , D h W t , s ) L2{RS+m + u£((blbt) t D h W, r , D h W t , s ) L2{R3+m 
+ u j\D h (m) t W* , D h W tlS ) L2(R 3 ;R3) + j\D h C trks W k , r ,D h Wi, s ) l2(r3 _ ;R) 
< C 7V(« , /) 2 + j (F u , D_ h D h W t ) L 2 {R 3 + . R3) + (D h H it , D h W ui ) L 2 {R z +m 

+ f {F 2t ,D_ h D h W t ) L 2 {R 3_ m + (D h H it ,D h W t ,i) L2{R3 _ m . (10.17) 





Since the estimation of the integrals with the indefinite sign in this inequality does not 
create any new difficulty with respect to the estimates that we have obtained in the previous 
subsection (they are even easier since the more difficult integrals 1$ and Iq do not have an 
analogue here), we provide the details in the appendix. With 5 > to be fixed later, this 
leads us to 

J o \\D h VW t \\ 2 L2{R3+ . R9) + sup \\D h VW\\ 2 L2{Km 

< [ CS(1 + C(M))N(u , ff + C5 1 (l + TC{M) + N(u , ff) ] \\(w, q)fz T 

+ Cs 1 N(u ,f) 2 + C Sl Ti \\(w,q)\\l+C s C(M)T$ (\\(w, q)\\ 2 + N(u , ff). 



We remark here that the estimates obtained in this section could have been performed with 
t € (0, T) generically replacing T; this explains the presence of 

sup || -DfcVW || ;R 9) 
[o,T] + 
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on the left-hand side of this inequality. As this inequality is independent of h, we then 
deduce that 



C llVoVWiH^ 8) + sup || Vo V^||2 2 

JO y + ' \0,T] 



[0,T] 

< [ C5(l + C(M))N(u , ff + Ctfi(l + TC(M) + N(u , ff) } \\(w, q)f Zr 

+ C Sl N(u ,ff + C Sl T 1 * \\(w,q)\\z T + C S C(M)T l 4 (\\(w, q)fz T + N(u , ff), 

and thus for the trace, where we will denote for notational convenience M. 2 = {x% = 0}, 

11^0^*11^0.5(152.^6) + SUP 1^0^11^0.5(^2.^6) 

<[C5(l + C(M))N(u ,f) 2 + C5 1 (l + TC(M)+N( Uo J) 2 ) ] \\(w,q)fz T 

+ Cs 1 N(u ,f) 2 + Cs 1 Ti \\(w,q)\\ 2 z T +C S C(M)T^ (\\(w,q)\\z T + N(u , ff), 

which implies that 

f T 2 2 

< [ CS(1 + C(M))N(u , f) 2 + CS 1 (1 + TC(M) + N(u , ff) ] \\(w, q)\\z T 

+ C 5l N(u , f) 2 + C Sl T\ \\(w,q)\\l T + C S C(M)TZ (\\(w , q)fz T + N(uq, f) 2 ). 

Since this has been done for any W = Q w o ^/j, we then deduce by the finite covering 
argument and the fact that each is of class H 3 that 



f 1 

I ||~||2 ii~i|2 

lFt|lHi.5(r o; ]R3) + sup IMlHi'5(r o; R 3 ) 

< [ C8(l + C(M))N(u , f) 2 + ChO- + TC{M) + N(u , f) 2 ) ] \\(w,q)f ZT 

+ Cs 1 N( Uo J) 2 + C Sl Ti \\(w,q)\\ 2 z T +C S C(M)Ti (\\(w, q)\\l T + N(uq, ff). (10.18) 

Elliptic regularity for the Stokes problem (see (for t £ [0, T] considered as fixed) 

-vA[w\ o fj- 1 ] + (q t o fj- 1 )* = -w\ t ° T 1 + Fl ° rf 1 + afw\ k ) t o f,' 1 

~ [2?>j Q[t V' 1 ~ [{ai)tq,k } o fj" 1 
+ u[(a J l af) t w' 1 , k ],j or)' 1 in fjfaSll) 
dw(w t o f)" 1 )^, ■) = -{(a{) t w\j } o Tp 1 in fj(t, fijj) 
w t o ff~ (t, ■) = on fj(t, Oil.) 
w t o fj~^ (t, ■) = w t o r) _1 (t, •) on fj(t, T ) , 
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then implies with the L°° in time estimate (j6.8|) of r\ (and thus of fj) into 

\\wt o ^~ 1 (*)lltf. W(t ,n/ ) ;R3 ) + 11% o rHtJIIfli^nJ)^) 

< C[ll " «4 V' 1 + F l t o f,' 1 - [(a*) t q, k ] o ff x + v[(aj afitw**],, ^W^^^f). 
+ life afw, k )t\\ 2 L2{mn ^ m + IK o 7? _1 (*)llifi'S(77(t,ro);K3) 

+ HK ! ^]*°r 1 n" 2ff(t)n / ) . R) ] • 

Thus, still with (EH) and (HT71) . 



II^(*)IIh2(^ ; k3) - C l|Vu> t (i)|| i4(n /. R9) + nwwiiffi^/. 
^ c I IKIIi^n^) + HVwII^,^/.^ + l|V™|| i4(n /. R9) + l|V^|| L4(n /. R9) 

+ II V 9Hl4(^;M3) + H^IIl 2 (0;R3) + |Klbi.5(r ; K 3) + Vf \ \ q t 1 1 L2 (n / ;R) ] 

from which we immediately infer that 

2 i / 11=; /+M|2 



11^)11^^3) + j o mw HHU f f 

T r-T 



<C[j \\wttf L2{n f. R3) + f ||/t|lia(n;R3) + / Rlllfi-5(ro;lR3) 
t/ ^ ** ** 

+ Ti Ml VxD! f Hl{nlm + T jf || Vw tt f H1{nlm + ^ || VwtWl 



+ Tl Nlv^olli 1( ^ ;R9) + t£ ||v^||^ 1(n , ;R9) + jT \\Vw\\ 2 H2(nf . 



+ T\[N(u jf + Tj Q \m\ 2 L2(<M) + 1 
Thus, with the trace estimate ()lU.18j) and (|1U.16|) 



r i-t 



2 r ~ 2 

H 2 (n^;R3) + y o Il^llj3-l(n^;][ 



< [C5(l + C(M))iV(no,/) 2 + ^i(l + rC(M) + iV( U o,/)') ] \\{wA)\\z T 

+ C Sl N(u J) 2 + CsC(M)T 1 * (\\(w,q)\\l T +N(u ,f) 2 ) + C Sl T^(w,q)\\z T ■ (10-19) 

Similarly, the classical elliptic regularity theory for the elasticity problem (for t £ [0, T]) 

-[c'^.i]^ = -w) tt + /t in fig 

w(t, •) = •) on T = 8Qq , 
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impl ies tha t ||^t> || J F/2(^g. M 3) < C [ \\ - w tt + ft \\i?(n*;M?) + IHI//i-5(r ;]R 3 ) ] > which with (|lU.18j) 
and (jlU.16f) provides us with the estimate 



SUP |H|#2 (n g J 

[0,T] 



< [ C5(l + C(M))N(u , ff + Cdt(l + TC(M) + N(u , f) 2 ) ] \\(w, q)\\z T 

+ C 8l N(u Q ,ff + CsC{M)T\ {\\{wA)\\l T + N{u J) 2 ) + Cs 1 T^{w,q)fz T • (10.20) 



10.3. Estimate of w independent of the regularization of a. 

Just as in the previous subsection, W again denotes W{ = C 2 w o ^j, where recall that 
^fi denotes the ith chart. 

Choosing <p = ^ -h^h^ -hPhW m the variational formulation ()9.64j) . we then find that 



±\\D_ h D h W(T)\\l 2{R3 . R3) + v j\blb s k D_ h D h W, r , D_ h D h W, s ) L2{Km 
+ l(C irks D_ h D h £ W k , r , D„ h D h £ W\ s ) L2(R3 _ ;R) 

- F (D„ h D h (bjQ), D-hDhW\j )l2( K 3 ;R ) 

>/ 

+ u [ T (D_ h D h (b r k b s k )W, r , D_ h D h W, s ) L2{R s +m 

J 

1 fT 

~^2 U (D ( -i)ph(b r kb s k )D(-i) Ph W, r , D_ h D h W, s ) i2(R 3 ;R3) 
+ j\D_ h D h [C* rks ] £ W k , r , D„ h D h W\ s ) i2(R3 _. R) 

- z2 / ( D (-i)ph[C irks ] D(_ 1 y h W k ,r ,D_ h D h W\ s ) L 2 (R 3 ;R) 
p=0 Jo Jo 

< C iV(n , /) 2 + / OD-ftFi, I^C-fclWWs ;R3) 
jo + 

+ / {D-h.DhHi, D- h D h W,i )i2( R 3_ ;R3 ) 

" 

+ / (D_ h D h F 2 , D_ h D h W) L 2 {R s_ m + (D_ h D h Ki, D h W,i ) i2(R 3 ;R8) . (10.21) 

»/ 
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Similarly as in the previous subsection, the estimates provided in the appendix yield (with 
5 > to be fixed later) 

\\D_ h D h VWf L 2 (M 3 + . K9) + sup \\D_ h D h V jf W\\ 2 L2(R3+ . R9) 

< C5 (l + C(M))N(u ,f) 2 \\(w,q)\\ 2 z T +CN(u ,f) 2 
+ C S C(M)T* (\\(w,q)\\z T +N(uojf) . 

As this inequality is independent of h, we deduce just as in the previous section that 



H*llH 2 -5(r ;iR 3 ) + s q u P II J o *ll/f 2 ' 5 (r o; R 3 ; 



< 



CS (1 + C(M))N(u , ff || (w, q)\\z T + CN(u , f) 



+ C S C(M)T* (\\(w,q)\\ 2 z T +N(u J) 2 ) . (10.22) 

Elliptic regularity for the Stokes problem (see |T2]) for t £ [0, T] 

-uA[w o ?? _1 ](t, •) + V^ofp 1 )^, •) = -wtofj- 1 +/ + ^,jOj ? " 1 (woj ? ~ 1 ), i 

- (a^-g) 77 _1 in fj(t,Ql) 
div(w o rf l )(t, •) = in fj(t, Q f ) 
(iiof'Kv) = on fj(t,dtt) 
(w o fi~ l )(t, ■) = (wo fi~ l ){t, •) on j/(t, T ) , 

then implies with (|6.8|) 

II* ° r 1 (*)ll H 3 W ( t ,n/) ; K3) + II? ° ^ l ^)Wmm,nl)m 

< C [ || -w) t o ? 7- 1 + / + ^, J or ? - 1 (u}o7 ? - 1 ), / || H1( . {t ^ );R 3 ) 

+ life ?) r^lffi^^n/j-R) + II* '7~ 1 (<)llH3-»«(t,ro);R 8 )] • 

Thus, with (ESI) 

H*(*)llH3(n/;M3) - C ll«K*)llwa.4(n£tt3) + ll?(t)llfla(n/; R ) - C Mit)\\ w iA { nf.«.) 

< C [ 11^11^1(0^3) + IHI W 2,4(fi/ ;R 3) + l|/l|i7l(0;R3) + 1 1 W \ \ H 2.5 (r o; R3)] , 
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from which we immediately infer, 

T r T 

2 i / II ~/j.M|2 







w Wllff3 (n /. K3) + y o iiywn H2(n /. 

T r-T pT 



< c[TN(u j) 2 + T 2 J \\w t \\ 2 H1{n f m +J o WfWhvtm+l W^H^m 



+ CTZ i\\ U o\\l Hn f.^ ) +Tj o ll^ll^ (fl /. R3) + j o NlH3 (n /. K3) ] 1 • 

Thus, with the trace estimate (jl0,22j) . 

T r-T 

II ~ II 2 4- / II ~ll 2 

<[C5(l + C(M))N(uoJ) 2 + C5 1 (l+TC(M) + N(u ,f) 2 ) ] \\(w,q)\\z T 

+ C 5l N(u J) 2 + CsC(M)T l * (\\(w,ml T +N(u J) 2 ) + C Sl T^\\(w,q)\\ 2 ZT . (10.23) 
Similarly, elliptic regularity for the elasticity problem (for t 6 [0, T]) 

_\Jjkl [ ~k 



: f w k ,i],j = -w t + f in Q s 
Jo 

[ w(t,-) = [ w(t r ) on r = dn s , 

jo Jo 



implies that 

*llH 3 (ng;R 3 ) <C[\\-W t + /||Hl(ng;R) + II / w l|# 2 - 5 (r ;R 3 ) ] > 

J 

which with (j!0.22|) and (|10,16|) provides the inequality 

sup || / u;||^ 3 (n«. R 3) 
[o,t] Jo v ' 

< [ C8(l + C(M))N(u , f) 2 + CStil + TC(M) + N(u , f) 2 ) ] \\(w,q)fz T 

+ C Sl N(u J) 2 + C s C(M)T^ (\\(w,ml T +N(u J) 2 ) + C Sl T^\\(w,q)\\ 2 ZT . (10.24) 

10.4. Existence and uniqueness of a smooth solution for the non-regularized 
system flTTTI). We now infer from (fTTHnT) . (fT0~T31l . (tTu~2771) . (HTmi) and (IT7H311 that 

ll(^)lll T 

< [C5(l + C(M) + iV( U o,/) 2 ) + C75 1 (l+TC'(M) + iV( UO) /)2 ) ] \\(w,q)\\l T 
+ C Sl N{u^f) 2 + C & C{M)T\ (||(u;,g)||| T +iV(no,/) 2 ) + C7 (5i r3||(*,g)||| T , 
this inequality being independent of the smoothing parameter of a. 
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We will call the constant C in this inequality C\ to indicate that at this stage it is a 
constant given by our successive estimates which, for the sake of conciseness, we have yet 
to make explicit. 

First, we fix b\ so that 

Cifc < ~ anddSxNiuoJ) 2 < J. 

o o 

The constant Cg x becomes thus determined, and we have that 

\\(w,q)fz T < [CxS (l + C(M))N(u J) 2 + 0^0(^1)] \\(w,q)fz T 
+ [C Sl + CsC(M)] T*\\(w,q)\\ 2 z T + C Sl N(u , ff 
+ CsC(M)N(u ,f) 2 T^ + ^\\(w,q)\\ 2 z T - 

Now let 

M = sup(M , 2 [d + C Sl ] N(u , f) 2 ) . (10.25) 
Consequently, C(M) becomes a fixed constant. Now, let us fix 5 > small enough so that 

\\(w,mz T <[l+C 1 d 1 TC(M) + [C 6l +C s C(M)]Tl] \\(w,q)\\z T 

+ C S C(M) T3 N(u ,f) 2 +C Sl N(u ,f) 2 + ±\\(w,q)\\z T ■ 
Now, let T S (0, Tw) be small enough so that 

|ll(^g)lll T < \ \\(w,q)fz T + Ci N(u ,f) 2 + C Sl N(u ,f) 2 , 

which implies 

\\(w,q)\\z T <M . (10.26) 

Henceforth, we revert to our original notation, denoting w and a by the sequential nota- 
tion w n and a n , respectively. The uniform bound (|10.26j) ensures the existence of a weakly 
convergent subsequence {w a r n \ , q a i n \ ) in the reflexive Hilbert space Yt such that 

(w a ( n ),Qa(n)) ~* (w,q) in Y T . 
The usual compactness arguments then provide the strong convergence 

K(n),<Z«x(n)) - K<Z) m L 2 (0,T;tf 2 (^;R 3 )) x L 2 (0, T; H 1 ^^; R)) . (10.27) 
Combined with the strong convergence 

o n o in L 2 (0,T;H 2 (Q f Q ;R 9 )) 

(which follows from the mollification process), the Sobolev embeddings provide the strong 
convergence 

a n alVw n -» aa T Vw in L 2 (0, T; L 2 (ft£; R 9 )) , 
a T n q n ^a T q in L 2 (0, T; L 2 (^; M 9 )) . 
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We then deduce from (|9~2Tj) that for each 4> E L 2 (0, T; FL\(S1^\ R 3 )), 
J (^»L 2 (n;R3) dt + vj (a r k w, r , a s k (j), s ) L2{n f. RS) dt 

+ J (C ijkl J ^ fc ,«,^,i)L2(^;K) dt + j Q (?> 4<^ ) L 2(n£;R) dt 

= / l ? !, ) L 2(n^R3 ) + (/, 0)L2(ng;R3) <ft , (10.28) 



which combined with the fact that, from (|1U.27[) . for all t E [0,T], w(t) E V v (t), proves that 
w is a weak solution of (|7.1|) . 

Since w E L 2 (0, T; Hq (Q; M 3 )) we infer the uniqueness of a solution in to this system 
in the same classical fashion as for the solution w of the regularized problem. 

Moreover, it is also immediate that we have from ()1U.26|) the estimate 

HItW < M . (10.29) 

This concludes the proof of Theorem 17. II 

Henceforth, M is given by (|1U.25|) and T is chosen such that (jlU.26f) holds. 

11. THE FIXED-POINT SCHEME FOR THE NONLINEAR PROBLEM 

We will make use of the Tychonoff Fixed-Point Theorem in our fixed point procedure 
(see, for example, jH]). Recall that this states that for a reflexive separable Banach space 
X, and C C X a closed, convex, bounded subset, if F : C — ► C is weakly sequentially 
continuous into X, then F has at least one fixed-point. 

With the quantities M and T being defined as in the previous section, we make the 
following 

Definition 11.1. We define a mapping Qt from Ct(M) into itself (from estimate \10.29\) ). 
which to a given element v E Wt associates w E Wt, the unique solution in Yt of \7.1\) . 

We next have the following weak sequential continuity result. 

Lemma 11.1. The mapping Qt associating w to v E Ct(M) is weakly sequentially contin- 
uous from Ct(M) into Ct(M) (endowed with the norm of Xt)- 

Proof. Let (i> p ) p gN be a given sequence of elements of Ct{M) weakly convergent (in Xt) 
toward a given element v E Ct(M) ( Ct(M) is sequentially weakly closed as a closed convex 
set) and let (u CT (p))pgN be any subsequence of this sequence. 

Since V^(T) is compactly embedded into L 2 (0, T; H 2 ^^; M 3 )), we deduce the following 
strong convergence results in L 2 (0, T; L 2 (£7q; R)) as p goes to oo: 



- 4 a l > ( 1Lla ) 

l(4U4) P }j^ (<W)j, (n-ib) 

(af) p ^at (11.1c) 

forip^fon. (ll.ld) 
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Now, let w p = ®t(v p ) and let q p be the associated pressure, so that (q p ) P £i>s is in a bounded 
set of Vj(T). Since Yt = Xt x Vj(T) is a reflexive Hilbert space, let (?*V(p), Qa'(p))peH be 
a subsequence weakly converging in Yp toward an element (w,q) £ Yt- Since Cr(M) is 
weakly closed in Xt, we also have w S Ct(M). 

We can then infer in a similar fashion as for the proof of Q1U.28|) in the previous section 
that for each <p £ L 2 (0, T; H^(n f Q ; " 



(w t , <£) L 2 (n . R 3) dt + v (a r k w, r , a s k (j), s ) i2(n /. R3) * 



+ / / w k ,i,(f) i ,j) L 2 {n s j . M) dt+ j (q, a{4> k ,i) L2{n f m dt 



(f°V^) L 2 (n f. R3) + {f, 0)£3(ng;R3) 



which combined with the fact that, from (jll.lj) . for all i £ [0>^1> 10 (^) £ shows that 

w is a weak solution of Q7.1JI in Cy(M), i.e. it? = 0r(u). 

Therefore, we deduce that the whole sequence (@T(v n )) n&N weakly converges in Ct{M) 
toward Qt(v), which concludes the proof. □ 

12. Proof of theorem Iq~T1 

The mapping being weakly continuous from the closed bounded convex set Ct(M) into 
itself from Lemma 111. II we infer from the Tychonoff fixed point theorem (see for instance 
[S]) that it admits (at least) one fixed point v = Q(v) in Ct- Moreover, since T < Tm, 
Lemma If). 51 ensures us that there is no collision between solids or between a solid and d£l. 
Thus, (v,q) is a solution of (|3.2j) . Note that the continuity of the Lagrangian velocities 
v* = v s at the interface To is ensured by our functional framework, since (y, q) G Xt 
implies v £ L 2 (0, T; Hq (0; M 3 )), which provides the equality vf = v s in H2 (r ; M 3 ). 

13. Uniqueness 

Uniqueness will be obtained under stronger assumptions than the ones used to establish 
existence, for reasons that will be explained hereafter. 
If (v,q) G Yt is another solution of (|3.2j) . then, 

(v - v)\ - v(aja k (v\ k -v\ k +a k (q, k -q, k ) = 5f in (0, T)xflJ, (13.1a) 

a k (v-v)\ k = 5a in (0,T)xft£, (13.1b) 

(„-«)*-[(«« / \v-v) k ,i],j = in (0,T)x^, (13.1c) 
j 

v(y* - vf )\ k afajNj - (q - q)^ = c^ kl [ (v s - v s ) k d Nj 

Jo 

+ 5g i on (0,T) x T , (13. Id) 

v-v = on (0,T)xdn, (13.1e) 

u-u = on {0} x n, (13.1f) 
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with 

6f = -v((aj4 - aj~4) v\ k ))„ +fo V -fofj + (-a* + S$)q, k in (0, T) X fi{ , (13.2a) 

5a = (of - 04)5*,* in (0,T) x fij , (13.2b) 

5g i = u(v f ,{ afajNj - v f ,{ afajNj) - q(aj - a{)Nj on (0, T) x T . (13.2c) 

If we view this problem with v — v as the unknown velocity and q — q as the associated 
pressure in the fluid, this problem looks similar to the linear problem (|7.1|) ; it is tempting 
to conclude that similar estimates as in the study of the regularity of (|7.1|) would yield 
a differential inequality that would provide uniqueness. It appears, however, that such a 
procedure fails because due to the Dirichlet boundary condition on d£l for the velocity, 
we are not able to get the necessary information on the second derivative of the pressure 
function. Such information is crucial since Sftt contains qu, which makes the second time 
differentiated problem impossible to estimate. 

For this reason, we will need to impose more regularity on the data and forcing functions, 
so that we, in turn, have enough information on q tt , which will then be viewed as a coefficient 
in the study of the regularity of 1)13. 1|) in Yt- 

We first update the functional framework. Let us denote 

Vf{T) = {u G L 2 (0,T; tf 4 (^;R 3 ))| u t G Vf{T)}, 
V S \T) ={uG L 2 (0,T;F 4 (^;R 3 ))| u t G V S 3 (T)} . 
Let us the denote the reflexive separable Hilbert space 

X% = {u G L 2 (0,T;fZ3(fi;R 3 ))| («/, / u s ) G Vf(T) x V S 4 (T)} , 

Jo 

endowed with its natural norm 

|2 _ ||„, 1 1 2 1 1 1 . 1 1 2 , 1 / „,l|2 



\m\x% - \\ u t\\x T + 11^11x2(0^4^/^3)) + II J Q u \\lHo,T;H^ s ;- 
In a similar fashion, we introduce 

Y? = {(u,p)\ P eL 2 (0,T;H 3 (n f ;R)), p t G L 2 (0, T; H 2 (Qq; 

' 



endowed with its natural norm 

\\(u,p)\\% =\n 2 X u + iipii^2 (0>Tifla( n/ iR )) + WptWl^o^H^uim) 

+ \\Ptt 1 1 L 2( 0tT . H l( Q f. R ))- 

We will also need 

W$ = {u G X%\ uut E L°°(0,T;L 2 (n;R 3 )), [ u G L°°(0, T; H 4 (n s ; ] 

Jo 

u G L°°(0,T;ff 3 (Og;K 3 ), u t G L°°(0, T; H 2 {tt s ; R 3 ) 
Uu E L°° (0,T; H 1 (Qq;M. 3 )} , 



<-;o 
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endowed with its natural norm 



u 



■ 2 iiii2 II ||2 11/ ||2 

lw™ = IMIx« + ll n «dlL°°(0,T;L 2 (C;IR 3 )) + II / n llL°°(0,T;_ff 4 (ng;R 3 )) 

J 

+ IMlL°°(0,T;.ff 3 (fig;]R 3 )) + II U t II L°° (0,T;_H" 2 (Qg;M 3 )) + II U « II L°° (C^T;// 1 (fig; 

as well as Z^(T) = {(u,p) G ijf| u G W^} endowed with its natural norm, 

IIKP)I|I| = NIw- + IIPllL»(o,T;H3 (n /^)) + H** llL8(0,T;fl3(n^;R)) 

+ IIP tt lli2( 0) T;ij-i(n/ ; R))- 

We can then define the convex set Cw-(T) in the same fashion as Ca/(T), with W-p 
replacing Wt and with the additional condition wu(0) = w-i where W2 has been defined in 
(jTHl) . 

We can then prove, in a way similar to the proof of Theorem 17. II (with the introduction 
of the penalized problems, time differentiated three times now) that the following holds: 

Theorem 13.1. With the same assumptions as in Theorem \5. 1\ and under the supplementary 
conditions Q of class H 4 , of class H 5 , the initial data uq G H 7 (Uq,M. 3 ) n H 4 (Qq;M. 3 ) n 
Hq(£1; M 3 ) n l? div j, /(0) 6 £f 5 (f2; 1R 3 ) ; satisfying the supplementary compatibility conditions 
( recall the assumption of Section \$) 

[(u[Vw f 2 Nf + 2uwl,{ (afaj) t (0)N J+ u4,l (a? af) tt (0)) 3 =1 ] tan 

= [{c^ kl w{, k l N 3 )U ] tan on r , 

w 2 G f#(fi;R 3 ), 

and the supplementary assumption on the forcing function that 

f G L 2 (0,T;# 3 (O;IR 3 )), f t G L 2 (0, T; H 2 (fl; M 3 )), / tt G L 2 (0, T; ff^fi; M 3 )) 

/ ttt GL 2 (0,f;L 2 (^;IR 3 )), 

we /iaue the existence of T > suc/i i/iaf i/iere exists a solution (v, q) G of \3. ty) . 
Furthermore, v G Cj^(T) for M appropriate. 

We can now get estimates for (|13.1|) which will give an appropriate differential inequality, 
in the space Zt used to prove Theorem 17. 11 We notice that this problem is similar to (|7,1|) 
except for the divergence-type condition which is not set to zero, and the boundary forcing 
on the interface. 

The Neumann forcing does not give any specific difficulty, and can be handled without 
modification of our previous estimates. 

The divergence- type condition does not bring any difficulty either because we do not need 
to establish the existence of a solution to ()13,1|) . since it comes de facto from the definition 
of v and v; we can directly use this condition in the steps where we obtained e- independent 
estimates for w#, w t and w. We also do not have to regularize the coefficients, since the 
regularity of w is a given. Those three steps would provide us in the same fashion as for the 
proof of Theorem 17.11 with the appropriate estimates to be made precise later. Note that 
this process works because the right-hand side of the divergence condition for w in (|13.1|) 
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has (roughly speaking) the term V77 — V77, which has one time derivative less than the 
right-hand side Vt> — Vt> (the term Vi; on the right-hand side being viewed as a coefficient 
whose regularity is given by Theorem 113. 1|) . 

We are now in a position to state the uniqueness result. 

Theorem 13.2. With the same assumptions as in Theorem, \1H.1\ and with the additional 
assumption that there exists K > such that 

Vt < f, V(x,y) eflxfi, 

\f(t,x) - f(t,y)\ + \Vf(t,x) - Vf(t,y)\ + \f t (t,x) - f t (t,y)\ 

+|V/ t (t,z) - Vf t (t,y)\ + \ftt(t,x) - f tt (t,y)\ + \V 2 f(t,x) - V 2 f(t,y)\ < K \x — y\, 

(13.3) 

i.e, f , Vf, V 2 f , ft, V ft and ftt are uniformly Lipschitz continuous in the spatial variable, 
then the solution is unique. 

Proof. With those assumptions on /, we have for the forcing f on— f of] appearing in 1)13. If) 
an estimate 

11/ 77 — /o v\\ L 2 {QtT . H 2 {n f. m + v)t\\ L 2 {0tT . HHn f. R3)) 

+ \\(f°V-f° v)tt\\ L 2 i0yT . L 2 {Q f m) < C \\v - ri\\vf{T) ■ 

The other terms associated to 5f, dg, 5a have the same effect in the integral estimates for 
wu, wt and w. This leads us to 

VtG(0,T), \\( v -v,q-q)\\z t < Ci \\v ~ v\\v*(t) ■ (13-4) 

where the constant C\ depends here on the same variables as the generic constant C as well 
as on the initial data. This thus implies 

Vt€(0,T), \\v-v\\ V 3 {t) < d Wri-fjWys^, 

from which we infer 

Vt€(0,T), ||w-i3|| V 8 (t) < Ci J \\v -v\\ V 3 (t) <dt ||u-«||v3( t ), 

which shows that for T\ = we have v = v on [0, T{\. We can then iterate this, starting 
from the initial time set at T\, which gives in a similar fashion, since v{T\) = v(Ti), that 
v = v on [Ti, 2T\]. By induction, we get v = v on [0, T]. □ 

14. CONCLUDING REMARKS 

Whereas the fluid-solid interaction is indeed a moving interface problem, it appears that 
the methods for its analysis differ drastically from the classical methods developed for the 
Navier-Stokes fluid interfaces independently by Solonnikov (see |2()| and references therein) 
and Beale j^]- 

First, our functional framework scales in a hyperbolic fashion in both the parabolic (fluid) 
and hyperbolic (solid) phases. 

Second, whereas the fixed-point problem (|7.1|) is inspired by the classical fixed-point 
problem used in parabolic-type interface problems, the Fourier transform technique used to 
get regularity in parabolic theories requires the introduction of the problem with constant 
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coefficients (for which one does get explicit solutions), with the forcing functions containing 
the difference (small in a neighborhood of a point on the interface) between the actual 
coefficient and this constant coefficient. Whereas this procedure is contractive for parabolic 
problems, the hyperbolic part is problematic in the sense that the difference between the 
actual and the constant hyperbolic viscosity is not regular enough to get these contractive 
estimates (those coefficients are not constant after the truncation and change of variables 
to the full space problem). 

Third, whereas energy methods without the use of Fourier techniques are indeed known 
for incompressible fluid interfaces (see for instance [8]), the highest-order time derivative of 
the pressure is known in L 2 (0, T; L 2 (f2; M 3 )) in that case, which allows the use of an iterative 
method from the constant-coefficient problem in energy spaces similar to the ones described 
in [20J - I n the fluid-solid problem, the knowledge of the highest order time derivative of 
the pressure is not known, which prevents such an iterative procedure from the constant- 
coefficient problem to get regularity. Instead, we are forced to work directly with the 
Lagrangian formulation (|7.1jl . which requires the introduction of the penalized problems 
for reasons explained previously about the pressure. In turn, working with the Lagrangian 
formulation (|7.1j) requires us to first smooth the coefficients, and then to obtain estimates 
independent of the smoothing parameter by using interpolation inequalities. 

Fourth, we clearly identify in our method the central and sufficient role of the trace of 
the velocity on the material interface Tq, whereas classical regularity results in interface 
problems involve the study of the regularity in the interior. 

Fifth, once again regarding the pressure estimate, obtaining a contractive fixed-point 
scheme does not seem possible for the hyperbolic-parabolic problem (even with data arbi- 
trarily smooth), whereas it is indeed the most well-known method for the parabolic interface 
case. Note, however, that this later point is associated to the incompressibility of the fluid 
and does not seem to appear without this constraint. 

This last remark is not without consequences for the numerical analysis of the problem, 
which we shall develop in future work. As for the question of the convergence of solutions of 
certain regularized models considered by other authors, it seems that the evidently natural 
approach of taking an elasticity law with a finite number N of modes introduced in 
and letting N —* oo leads to some difficulties, as there is no elliptic operator for the discrete 
elasticity problem for which one may use H 3 regularity independently of the number of 
modes. On the other hand, it can be shown that the addition of a hyperviscosity to the 
solid problem (similar in spirit to the hyperviscous plate problem introduced in |S]) would 
indeed converge to the solution of the actual problem as the hyperviscosity parameter tends 
to zero, since we can apply the methods constructed here to this family of problems and 
obtain estimates that are independent of the hyperviscous parameter in the correct norms. 

Appendix A. Some additional estimates 

A.l. Estimates for (|10.17|) . In this section, 5 > is assumed given and we now proceed 
to the estimate of the terms of Q10.17J1 whose sign is indefinite. Recall that from (jfi.^j) . a, 
at, and thus b, b% are controlled respectively in L°°(H 2 ) and L°°(i? 1 ), independently of the 
regularizing parameter n associated to a. 
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Step 1. Let Ji = / (D h b( Q t , D h Wij ) i2(R 3 + . R) . Then, 

J 



f T ~ 2 [ T 5 15 

~Jo H jD/lVTy 'l^ 2 K;iR 9 ) + c s c ( M ) J \^L 2 (M 3 + ;M)\\Qt\\HHn 3 +] I 
< 5 £ \\D h VW t \\l 2{Rlm 

+ C S C(M) T\ [sup ||Q t ||£ 2(R3 . R) + f 

CO.T) V +' ' JO 



(0,T) 

From (|1U.1|) and the definitions of and Q, we then infer 
fT 



T 



2 



|Ji| < f 
Jo 

+ C S C(M)T* [N(u ,f) 2 +T^ llftll^ (Q / ;R) + sup|l^lli 2 (n;R3) 



[0,21 

/•'/ 



JO 

+ sup ihIhi^rs) + y o ii^n ffl(n /. R) ] • (a.i) 

Step 2. Let J 2 = / (AAi Q, D hWl,j ) L ^m?, ; r)- Similarly as for J x , 

|Ja|<* y HA»VWt||£a(tt3 ^9) + Cfc y \\D h b t \\l 2(R3+ . R9) \\Q\\^ 1A{R ^. R) 
< s£ \\D h VW t \\l 2{R3+iR9) + CsC(M) jf ||Q||^ (R 3 +;R) ||Q||], 5 2(R 3 +;R) 

<^/ T H- t H 2 H 2( n/;R3) 

+ C S C(M) rJ [ i>T(«o,/)» +r jT llftll^^j + jT llffll^tp/.^] • (A.2) 
Step 3. Let J 3 = / (6 t - D h Q, D h Wt,j) L 2 (R 3 + . R) . Then, 

« 

\Js\<5 F \\D h VW t f Lm s R9) + C s f T \\btf LHR s m \\D h C» 2 
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Thus, similarly as for (|A.2|) . 

(■T 



2 

/.it 



|Js| < C<5 / 

+ C 5 C{M) T^iV(n ,^ + T^ T ||^Ji 1(Q , ;R) + ^ T ||,1^ (Qo/;R) ] . (A.3) 
Step 4. Let J4 = / (6] D^Qt, ^hWl 13 )l 2 (r 3 -r)- This term will require more care. We 





first notice that 



J 4 =/ p h Qt, D^Wlj})^.^- [ (D h Q t , D h V i W{, j (- + h)) LHRl]i 
Jo Jo 



'0 -/ 

which with the divergence relation (j9.66f) leads us to 



Ja = / (AiQt, ^>hdt)L2( R 3 ;R) - / (D h Q t , D h [b t lW\j]) L 2r R 3 m 
Jo Jo + 

- / (Z> h Q t , (• + /i) W ;R) . (A.4) 

Jo + 

f T 

For the first integral of this identity, j\ = I (DhQt, D^a^^ma -R) > we h ave 

J 

1 41 < &£ \\ D hQtf L 2 (Km + C s £ \\D h bA\ 2 r,^ .^\\w" 2 



^ W tllL2(R3.;R9)ll" ; ll W /i,4 {n / ;R 3 ) 
T ,-T 

+ C S J q l|- D ^Hi4 (K 3_. R 9 ) ||^tlli4 (n /. R 3 ) + C 5 ^ IHIwl.4(R3 + ;R9)ll V ^lli2 (0 / ;R9) 

+ C s ll 6 *ll^ (R 3. R9) ||Vu;||2 4( ^ ;K9) 



< s£ \\D h Qt\\\ HKm + C 5 C(M) T [sup \\M\ 2 HH<<m + sup 



T 







+ CsC(M) T [N(u ,ff + T £ \\w t f H2{nfom+ T£ \\w tt \\ 2 Hl{n f m ] . (A.B) 
/•T 

Next, for jf = / (D h Q t , D h \b t \W l 
Jo 

\J%\<5 \\D h Q t \\ 2 L2{R s m + C s I ||-D/ l ^|li2 (R 3. R 9 ) ||VPF||^ 1 , 4(K 3. R9) 



f T ~ 2 2 
+ C S ll^ll|4 (R 3. ]R 9)||-D? l VW|||4( ]R 3. R 9) • 
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Therefore, 

\Jl\ < $ \\D h Q t \\ 2 L2(R3+ . R) 

+ C 5 C{M) T\ [sup||V^||^ 1(R 3 ;R9) + ^ T ||V^||^ (R 3 ;R9) ] , 
which with the definition of W and Q provides 

\Ji\<cs^\m\i 1(nlm 

+ C 5 C{M)T\ [^o,/) 2 + T^ T ||V*J^ 1(Qo/;R9) + ^ T ||V*||^ K;R9) ] . (A.6) 

f T 

Similarly, for jf = / (D h Q t , D h b? W t \j (• + h)) L2(R 3 +m , 

J 

3 f T 2 f T ~ 2 2 

\JI\ < ^ / \\ D hQt\\ 2 L 2^. R) + C S / ||-D h 6||^ 4 ( R 3. M 9)||VW t || L 4( M 3. R 9) 

•/ « 

(A.7) 

Step 5. Let J 5 = f (D h (b r k b s k )W t , r (• + h), D h W t , s ) L 2 (R 3 + . R3) . Then, 



rp rp 

\M< sJ q \\D h vw t \\ 2 L2{R3+ . R9) + c s J o ll^(^ T )lli4 (R 3 ;R9) ||V^||i 4(R 3 +;R9) 

+ C S C{M) T\ [^oJ) 2 +T^ T ||V*dli 2( ^ ;R9) + ^l|V^|li 1( ,/ ;R9) ] • 

(A.8) 

Step 6. Let J 6 = / (D h (b r k b s k ) t W, r {■ + h), D h Wt, s )^ R 3_. R 3). Similarly, 

J 

|J 6 |< 5 f T \\D h VW t f L2{K . R9) + C s / T ||^(6F) t || 2 :2(R 3. R9) ||VW|| 2 yl ,4 (M 3 + . R9) 

+ C s C(M)rJ [Ar(»o,/) 2 +r^ T ||V»i < || 2 ff , (n , iItJ) +^ T ||V»f H2(n , iIt9) ] . (A.9) 
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Step 7. For J 7 = f ((b r k b s k ) t D h W, r , D h W t , s )L2 (ffi 3_. K 3), we have 

J 

+ C 5 C(M) T* [^oJ) 2 + T^||V*J^ K;R9) + ^ T ||V*|| 2 H2K;R9) ] . (A.10) 

For the next step, we introduce b\ > 0, which is different from S. 
Step 8. Let J 8 = / (Z^C^H^V , D ft W?„ ) L2(M 3_ 



An integration by parts in time gives 



h = - f (D h C irks Wj°, r , D h W i , s ) L 2( K 3 ;K ) 
Jo 

+ [ ( D h C irks W\ r (•), ^^(-), s ) L2(K 3_ m ft . 



Since fig is of class H 4 , 



and thus, 



\M<CT[ sup || VW t ||2 a(R3 _ ;R9) + sup ||W||^ (R 3 ;R3) ] + CN(u , ff 

+ C Sl [ \\VW(T) - W(0)Hi 2(M 3_ ;K9) + ||V^(0)||| 2(R3 _ ;R9) ] 
+ 5! sup H^WII 2 3 . R9) , 

[0,T] V - ; 



| J 8 | < CT [ sup ||Vu> t ||2 2( 9) + sup ||*||^ 2(Q » RS) ] + C Sl N(u , ff 

[0,T] ' [0,T] y ' 



+ C Sl T sup ||V«) t || La(ns R9) + C6i sup ||«)|| Ha(n .. R s ) • (A.11) 

[0,T] y ; [0,T] ; 
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Step 9. Let J 9 = f (F u , D_ ft D ft W t ) L 2( K 3 ;)a 8) . Then 

J 

\h\ < 6 [ T \\Wt\\ 2 HHK , R3) + C 5 / T ||(66 T ) t |li 4(K 3 ;K9) ||VTy||i 4(K 3 



+ C 5 J [K^II^MfRa^llVWill^a^j+CiiVCuo,/) 2 

f T ~ 2 2 /" T " 2 2 

+ <?5 ^ ll^lli 4 (M 3 + ;M 9 )ll«lli4 (c / ;M) + C5 ^ ll b lli- (R 3_. ffi 9 ) ||%||^ 2(n /. I 







+ C 5 C(M) T [iV( U0 , / )2 + T ^ T ||^||^ ( ^ ;R3)+T jj\\w u f 



Thus, with (tlTHl) . 

' l| ^ 2 (^o; 



|J 9 | < C6j T \\w t f u2 , nf ^+C s N(u J) 2 



+ C$C{M) T [^o,/) 2 + T^ T ||^||^ K;R3)+ r [\\wuf m{< ^ ] 
+ CgC(M) T[T jT llftll^jn/.^ + sup ||t%||i a( n ; RS) + sup INI^.rS) ]. (A.12) 
Step 10. For Ji = / (D h H H , D h W t ,i )l2 (R 3 . r3 ), we have 



|Jl0|<*/ IM^ (R 3. R3) +Q/ ll(^ T ) t |li4 (R 3 ;K9) ||V^||^ 4(R 3 ;R9) 



JO 



o 



+ C s [ fbb 
Jo 



T 



< C6 



T 



r 1 

+ C S C(M) T [N(u ,f) 2 + T 

Jo 



(A.13) 
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Step 11. Let Jn = J (D h F 2t , Df l W t ) L 2g S 3_.^ + J (D^ h D h K it , W t ,i ) L 2 (R 3_ ;R3 ) . Then, 
| J u | <CT [sup\\w t \\ 2 Hl(n s R3) + sup\\w\\%2 (n s m ] + CN(u J) 2 . (A.14) 

[0,T] °' [0,T] ; 

A. 2. Estimates for (|1U.21|) . 

As in the previous section, recall that from Q6.3JI . a, and thus 6, is bounded in L°°(H 2 ) 
independently of the parameter n associated to a. 

Step 1. For K x = / (£>_ A Z> fc (ft?) Q, D- h D h W\j . R)j we have 

■> o 



2 , /-» / II r> r> tll2 



^ J / T |^-^ W II^;k 9 ) + ^ C ( M )/ T |I Q I^Hk 3 + « I|Q|1 



T 

2 

tf 3 (£^ ;R 3 ) 

T /-T 



< / IN,,:.,,,' 



+ W(M) TUN(u ,ff + Tj o \\ qt f Hl{nfom + 1 IMI^ (n , ;R) ] • (A.15) 

1 ,T 

Step 2. Let K 2 = Y, ( D (-i)pfA D ( _ 1)Ph Q, D_ h D h W l tj Then, 
p=0 ^° 

1^1 < *jf P-^W||£ 2(M a ;R9) + C S £ \\D h bf LHR3+ m \\D h Qf LHKm 

- C6 l ^(olm 

+ CsC(M) T\ [^(no,/) 2 + r^ T || % ||^ lK;R) + ^ T ||< 2(no/; J . (A.16) 
Step 3. Let K 3 = f (6? D^ h D h Q, D„ h D h W\j ) L 2 {R s + ;R) . 

•/ 

;e that 

^3= [ T (D_ h D h Q, D_ h D h [^\ J ]) L2{Km 
Jo + 

+ ( D ~h D hQ, £(-l)Pfc&i ^-lJJ-fcWSj ) L 2 (R 3_. R) 

- f T (D^ h D h Q, D h D^ h [bl] W i , j ) L2(R3+;R) , 

•/ 



We first notice that 



MOTION OF ELASTIC SOLIDS INSIDE OF A VISCOUS FLUID 



which with the divergence relation (|9.66f) leads us to 

i-T 



K 3 = J (D- h D h Q, D_ h D h a 

1 t 

+ Y\ I (D- h D h Q, D { _ 1)Ph b{ D ( _ 1)Ph W 



p=0 • 

T (D_ h D h Q, D h D_ h [V>] W\ 3 ) LHK 



o 



We then have that 



r-T /-T 
l^sl < S I \\D- h D h Q\\ 2 L2(K 3 :K ,+C s IH| 2 ff2 m 3 :R 9j|*||L, /. 



o 

r-T 

+ C 5 ' 







+q y o T iHi^ ( R3 +;K9) iivwii 2 WM(M 3 

< C<5 ^ ||g||^ 2(n /. R) + C S C(M) T [ N(u , ff + T jf 

+ C 5 C(M) T\ [N(u ,ff +T j o T \\w t \\l 2{ ^ R3) + £ \\w\? H3{Qim ] • (A.17) 
Step 4. Let K 4 = ! (D_ h D h (bib s k )W, r , D_ h D h W, s ) L 2 { ^ 3) . Then, 

J 



fT 

r-T r-T 



T i-T 

T\\2 



\K A \< \\D_ h D h VW\\i 2{M3+ . w9) +C 5 j o ||6^||^ (M a ;R9) ||VW||^ M(R3+ 



T 



+ C S C(M) T\ ^(noJ) 2 + T^ T ||V^||^ 1(c , K9)+ ^ T ||V*||^ (n , K9) ] . 

(A.18) 



Step 5. For K 5 = Y] I {D^ iyh {blb s k )D^ 1)vh W, r , D„ h D h W, t 



t 

~ 1 1 2 



<^ / \\w\\..„~f- a 



i?3(fij$;R3) 

T i-T 



+ C s C(M)tUn(u J)* + tJ o \\m\\ 2 H2{n i m + 1 M 2 H3{n f m } ■ (A-19) 
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Step 6. For K 6 = j\n^ h B h [C irks ) J W k , r , D_ h D h W\ s ) L2(R 3_ ;R) , an integration by 



parts in time gives 

f-T 



K 6 = - f {D_ h D h [C irks ] W k , r , D_ h D h [ W\ t 
Jo Jo 

+ ( D- h D h [C irks ] f T W k , r , D_ h D h [ T W 
Jo Jo 



from which we infer from the H 4 regularity of 



o- 



K 6 <C I ||VW|| i 4 (M a ;R » ) p_ h Z) h y VW|| L2(K 3_. R9) 

hPh I VW|| L 2( R 3 ;R 9) , 

Jo 

leading us to 

l-Kel < CT [ sup \\wf H2 , Qam + sup || / wf HHna R8) ] . (A.20) 

[0,T] [0,T] io 

Step 7. Let 

K 7 = J2 [ T (D(-iyh[C irks ] I D { _ 1)Ph W k , r , D_ h D h W\ s ) L2(R3 . R) . 

p=Q Jo Jo 

An integration by parts in time gives 

K 7 = ~z2 I { D (-i)Ph[C ir s ] Dr_^ P +i h W k , r , D_ h D h l W\ s )l2(M? ;M) 
p=Q Jo Jo 

+ ( D ( _ 1)Ph [C irha ] D ( _ 1)Ph [ T W k , r , D_ h D h [ T W\ 

Jo Jo 



and thus from the H regularity of f2g, 

Ki< C £ ' \\D h VW\\ LHM 3_. R9) \\D„ h D h J V^|| i2(R3 _ ;R9) 

+ \\°hj ™\\v<R*_ja*)\\D-hD h J V^|| L2(R 3_ ;R9) • 

Therefore, 

\K 7 \ < CT [ sup \\w\\ 2 H2(n s R3) + sup || / w\\ 2 H s(mm 1 • ( A - 21 ) 

[0,T] [0,T] Jo 

Remark 21. T/ie -ff 4 regularity of fig is used only for proving \A . 20\) and \A.21\) . As a 

matter of fact, W ?J ' P for p > 3 would have been sufficient. 
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Step 8. Let K % = ! (D- h F u D h D_ h D h W) L 2 {R 3 + . R3) . Then 

J 

\K 8 \<5J^ \\W\\ 2 H3{R3+m + C 5 j ll^(^ T )lll4 {R 3. R9) ||V*||2 4( ^ ;R9) 
+ C s J q ll^ T |lloc (M 3. K9) ||^Vw)||^ (n ,. M9) + C N(u ,f) 2 

f T ~ 2 2 f T ~ 2 2 

+ Cs Jo " V6 '^ 4 ( R +; R27 )" 9 "i 4 (^« + Cs J ll 6 lli-(ffi^K9)ll v 9ll L2(Q /. M 3 ) 

+ C 5 C(M) T [iV(n , /) 2 + T ^ |kt||^ 1(Q ,. ffi) ] + C N(u , ff . (A.22) 
Step 9. For K 9 = f {D_ h D h H h D_ h D h W yi ) L 2 (R 3 ;R3) , we have 

J 

\ K ^\< S l \\ W \\ 2 H*(Rl;RZ)+ C sJ o ll^(^ T )lli4 (M 3 +;R9) ||V*| 

+ CsJ \\D_ h D h (bb T )\\ 



- -, - , ^1^1,4(^3) 

T 





+ C S I ll^ T |lL-(R3 ;R 9 ) |k||^ 2(nj f. M3) 



<S[ IH^ (n , ;R3) + C s C(M)T[N(u jf + Tj o |K|I^ ;K3) ] ■ (A.23) 
Step 10. Let 

K w = f (D_ h D h F 2 , D_ h D h W) L 2 {R3 _ m 

J 

+ / {D-hDhD_ h Ki, D h W,i ) L 2q SL 3_ ;K 3) . 

J 



Then 



K w \ < CT [ sup |H|| 2( ™ ;K 3) + sup || / w\\ 2 H z {n s R3) } + C N(u , ff . (A.24) 

[0,T] [0,T] JO 
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